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INTRODUCTION 

This volume describes the incremental computer from the operational point 
of view. In brief, the computer accepts real-time analog inputs and continu- 
ously modifies pertinent functions of the variables so as to yield updated out- 
put functions. For example, the fire control problem has as inputs the vector 
position and vector velocity of an attacker and computes as outputs the azimuth 
and elevation lead angles. Only the change of the output function is computed, 
and this is added to the initial function value to give the updated value. The 
computer operates by performing one basic computation (whose formulation is the 
basic algorism discussed later) about 64 times every drum revolution. The 
quantities necessary for a computation are read off the drum serially and out 
of the random access memory during a minor cycle. 

The detailed analysis of the mathematical operation of the computer and 
the points considered in programming and coding follow as outlined below. 

A knowledge of Section I is not strictly necessary to the practical opera- 
tion of the computer except perhaps for the section on errors. It is invaluable 
for a mathematical understanding of the computer action, however, and the spe- 
cial techniques discussed indicate interesting but less obvious possibilities 
of the computer. 
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1. MATHEMATICAL ANALYSIS OF OPERATIONS 

a. INTRODUCTION. - The application of digital techniques to the methods 
of solving mathematical problems normally associated with analog computers leads 
to a hybrid computer that has some of the advantages and some of the limitations 
of both the analog computer and the general purpose digital computer. The ERA 
incremental computer, as presently designed, carries out numerical computations 
by acting on the small, discrete changes of the input variables to compute the 
small, discrete changes in the output variables. In the memory of the computer 
are stored the orders for a number of individual steps (called "minor cycles") 
of numerical computation. During one major cycle, all the variables in each 
minor cycle undergo changes of one increment, i.e., plus or minus one. At the 
end of each minor cycle, the output of that cycle is stored in the random ac- 
cess memory as a plus or minus one. This stored output can be used to modify 
any of the variables in any of the several steps including the step from which 
the output is obtained. In any minor cycle, those numbers whose changes are 
determined by other minor cycles (through the random access memory) or by an 
input quantity are called independent variables, and those that are changed by 
the step itself are called dependent variables. 

An incremental computer is useful primarily because it is capable of doing 
fairly complicated operations within a small computer and it produces an up-to- 
date answer output once every major cycle, i.e., about once every 0.005 seconds. 

The incremental computer can be made small because its arithmetic section 
only adds, subtracts, and complements. These can be combined to form the more 

^ complicated operations within a single step, 

» 

w In principle, the incremental computer differs from a general purpose com- 

e ■ 

puter in that instead of performing a complete numerical operation at each step 
of each major cycle, it computes the change in that step's solution due to the 
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change in variables since the last major cycle. For example, instead of form- 
ing the sums U Q n _x + V Q T n-1 , U Q Q n + V T n , U Q n+1 + V T n+1 , etc., directly in 
successive major cycles, the incremental computer forms 

" O 0n + Vn = D oQn-l + v o T n-l + °»W. - «n-l> + V T n - Vl> 
"A[U Q„ ♦ V„tJ =0 o AQ n + V AT n , 

and A[D Q n+1+ V T n+1 ] = U„A Q n+1 + V AT fl+1 , etc. 

By restricting the change of each variable to a +1 or -1 each major cycle, 
these operations are accomplished by addition or subtraction. Through the use 
of these simple operations, appropriately controlled, the incremental computer 
can form sums, differences, products, quotients, square roots, integrals, deri- 
vatives, integrals with an inverse integrand, logarithms, and exponentials. 
Any one of these arithmetic operations can be done in a single minor cycle. 

It is basic to the use of an incremental computer that the error of one in- 
crement in any variable is small enough so that to neglect it would not cause 
serious errors in the answers. This implies that in five milliseconds the 
answers themselves change by a "negligible amount". This, in turn, implies 
that five milliseconds is a "negligible" period of time. In other words the 
incremental computer is useful for real time control of systems whose time con- 
stants are of the order of seconds or longer. In systems with relatively long 
time constants, the answers can be changed by large numbers of increments dur- 
ing a period of time short compared to the time constants. With longer time 
constants it is possible to represent the variables more accurately within the 



$ machine by making one increment a smaller part of the variable. The only prac- 
04 tical limit on the accuracy of the numerical solutions of the mathematical 
|n problem: programmed in the incremental computer is the length of time to effect 
the changes in the variables. 



In direct analogy to analog computers, the incremental computer can be used 
for the solution of implicit equations. The pairs of operations multiply and 
divide, integrate and differentiate, squaring and extracting a square root, and 
generating the logarithm and the exponential, can be considered as explicit- 
implicit pairs in the incremental computer. Simultaneous and differential equa. 
lions are solved by an implicit arrangement just as are similar problems in an 
ilog computer. They are fundamentally different, however, because the incre- 
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mental computer handles numbers digitally with the associated advantages 
numerically and electrically, 

b, MACHINE ALGORISM. - The overall problem to be solved by the incremental 
computer must be divided into a sequence of arithmetic operations or steps. 
The steps are coupled by the increments stored in the random access memory; the 
output increments from one step are the input increments to other steps. In 
the presently designed equipment, the heart of the storage system is the mag- 
netic drum. The commands and numbers required for each of the steps of the 
mathematical problem are stored sequentially on several tracks on the drum. 
Each step is acted on by the arithmetic section of the computer once each drum 
revolution, that is, one drum revolution corresponds to one major cycle. 

Let us consider the operations performed by the arithmetic section of the 
computer during a single step. During each step four binary numbers are read 
from the magnetic drum by four reading heads in the arithmetic section of the 
computer. These numbers will be denoted as Uj ^, V. , f Rj_i, and S, Three of 
these numbers, U j _, . v i_li an( * R i-1» are processed in accordance with the com- 
mands from the control section of the computer and the new numbers are returned 
to replace the previously read numbers in their respective places. The number 
S is a constant scale factor which remains on the drum. 
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At each step the arithmetic section forms the incremtal sums 

R i = R i _ l + UiAQj + V^jATi + SAPj - SAW| CI) 

To simplify discussion, the increment will be considered unity and all machine 
numbers will be considered integers. That is, the increments ^T^ 9 AQj, A Pj 9 
and A Wj, and the incremental changes in the variables Uj and Vj, namely AUj 
and AVj, can have the values +1 or -1, These incremental values are read from 
the random access memory from positions given by the commands interpreted by 
the control section. In many cases some of the increments are programmed to be 
zero for all major cycles of a particular step, i.e., they are not read from 
the random access memory. 

By summing a single step over n major cycles, the values of the variables 
during the nth major cycle can be determined. That is, 



Z CR| - R UI ) = I UiAQi + £ V^ATi + S £ A P, - S £ AW, 
i=l i=l i=l i=l i=l 

or R B - R = Z OjAQj ♦ I V^jATj + SCP M - P ) - SCW n - W Q ) (2) 



In the computer R Q = in every case. In all but special cases the incremental 
changes in the dependent variable are chosen to cause R n to go toward zero. 
When R n is sufficiently close to zero to be considered negligible, the step is 
"settled". When it has settled, the equation for the machine algorism for that 
is approximately? 

O = 1^ UjAQj ♦ I . Vj jATj ♦ S(P„ - P > - S(W n - W ) (3) 

In each case the increments are considered small. If we consider them to be 
differentials, equation 3 suggests? 



=.....'/ UdQ + / VdT ♦ S[p n - P ] - s[w„ - W ] (4) 






n 

/ 

o o 

c. QUALITATIVE DESCRIPTION OF ALGORISMS. - In order to determine exactly 
the course of an arithmetic operation performed by a step of the incremental 
computer, equation 1 must be used with the mathematics of finite differences. 
It is useful to use equation 4 in a qualitative way to approximate the answers 
actually obtained in the incremental computer. In applications, the situation 
is reversed; the problem is to find an incremental program which will approxi- 
mate the desired mathematics. It is hoped that the qualitative use of equation 
4 will clarify the principles of the arithmetic algorisms. 

Sum and differences; Set U and V constant and solve for W, Then from 
equation 4 

,n a n 

I 

o 



O = / o dQ ♦ / V dT + S(P„ - P ) - S(W„ - W ) 



O 



"n = s * P n - P + W 

Let the initial conditions satisfy 

„ _VojiMo + p 

W ~ S 

Then 

„ "oV Vn n 

w n = s + p n (5) 

That is, through the use of the sum algorism up to three variables may be added 
in the form 

X = Xj * 89X9 + 83X3 

where 83 and 83 must be rational numbers, &2 an< * a 3 are tne quotient of the 
positive or negative integers U and V Q divided by the integer S. 



ct and quotient: Set dT = dU, dQ = dV, a.ftd sojve- for W, Then 
n 

I 

o~ o 



J UdV ♦ / VdU + S(P n - p ) - s<w n - w > 
o J (UdV + 



»■ = w o ♦ P„ " Po + ~ ~ U s " VdU 

- w ♦ P P + U n V n " u o v o = w . p + HoZfi + p + Z& 
~ w o + r n " r o s ° ° S n s 

Let W = -p ♦ P ^ 

I 
V u 
Then W n = -J3U1 + P fit multiplication (6) 

o 

S(W n - P n ) 

or U n = , division (7) 

v n 

Clearly, from equations 6 and 7, these algorisms should more properly be 
called "Product and Sum" and "Sum and Quotient*' in general. 

Square Root: The square root of a variable can be derived from the incre- 
mental algorism by letting U = V - T =. <$. Then equation 4 becomes: 



n n 
= J VdV + J VdV + S(P n - P ) - S(W n - w ) 



G = V 2 - V Q 2 - S(W n - P n - W * P ) 



v 2 = S(W - P) + v 2 , S(W - p ) 



Let V 2 = S(W - P ) 

V 2 = S(W -P) ' 
and V = 7S(W - P) (8) 

(1) INTEGRAL. - In equation 4 the availability of integration is obvious, 
However, in a step doing a single integration it is sometimes advantageous to 
set U = V, X n = X = 0, and A W = AT. For more specific details refer to the 
quantitative description of integration. 



(2) TIME DERIVATIVE. - Differentiation with respect to time can be 
treated as the inverse of integration. However, , special care must be taken to 
make it stable. Differentiation is accomplished by letting the U variable be 
constant, the T variable be time, and dV = dQ and solve for V, Equation 4 be- 
comes: 

n „ n 
U dV + t 
o o 



o = /" u dv + / Vdt + s p n - p .. w n .♦ w 4 



Differentiating: 

- U o dt * V + S dt " ^t 

V = S^ (N -P) -U<$ (9) 

If U were zero, V would be proportional to the derivative of the function W - P, 
This algorism would properly be called "Addition and Differentiation", However, 
the last term with U ^ is necessary to obtain a stable differentiation algo- 
rism. This stability term can be compensated for. The accuracy of the compen- 
sation is determined by the number of minor cycles used. 

Integral with- a Reciprocal Integrand . 

Let U = V and solve for Q + T in equation 4: 

n M n 
VdQ ' + 
o o' 
. ^ / 

Differentiating: 

= Vd(Q + T) ♦ Sd(P - W) 
d(Q+T) s Sd(W^ P) 

Integrating: 

+ T = s/ <*^- P > (10) 

Usually Q = T and P = 0. Then equation 10 becomes 

Q = |/f (U) 



= / '■ VdQ + J VdT + S(P n - P ) - S(W n - W ) 



(3) LOGARITHM. - Reciprocal integration can be used to generate the 
logarithm. If we let V = W, equation 11 becomes; 

or Q = § log e W ♦ C (12) 

(4) EXPONENTIAL . - The exponential can be obtained from the logarithm 
by solving for W, Then 

|p *C 9 =log e W 

e S = e log e W 

W = e S ' ( (13) 

The foregoing has been a qualitative treatment of the. use of the KRU incre- 
mental computer. The following discussion is a more formal presentation in- 
cluding the specific restrictive formulas for the step and the choice of the 
sign of the increment of the dependent variable to make the value of R n go toward 
zero s 

d e QUANTITATIVE DESCRIPTION OF INCREMENTAL ARITHMETIC OPERATIONS, - The 
algorisms previously qualitatively described are specifically described by the 
choices of the increments and initial values. In any step there are six incre- ' 
ments to be specified (AU| 9 ATj, AVj, AQj, APj, and AWj), two initial 
values CU and V ) t and the scale factor (S). The initial value of Rj e that 
is, R Qt is always fixed at zero/ The correct combination of the programmed 
initial values, the programmed selection of the six. Increments In each step, 
and the wired-in machine algorism represented by equation 1 cause the incre- 
mental computer to perform the arithmetic operations described above as well 
as some special operations described later e Those variables of a given step 
that have their increments determined by other than the step itself, e 4 g eo by 
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other steps, are called independent variables. The increment of the dependent 

variable depends on the sign of the remainder, R n , of the step in question 

(which may vary from major cycle to major cycle) and usually on the non-varying 

sign of other variables or constants as well. For convenience of notation the 

signum function will be used. It is defined bys 

Signum xssgn xs + 1 for x > 

s - 1 for x < 

The programmer is free to choose the increments AU|, AVV APj,, AW, C 

AQj,- and AT| from any position in the random access memory. Each increment 

stored in the random access memory is either sgn Rj, fori one of the steps or u-.. 

s 9 n CV| n - Vj)- for one of the input quantities, 

(1) SUM AND DIFFERENCE . - To perform the previously mentioned sums 
M PoQn + V T n . 

set AUj = AV| = 

m PoQo + Vo o 

w o = """"s + P o €15) 

and AW i+1 = sgn Rj sgn S (16) 

The sense of equation 16 is such that the increments of W tend to reduce the 

magnitude of Rj, Equation I becomes : 

R i - R i~l = u o A Qi * V o AT i * SAP i ~ SAW i 
By summing this equation for n major cycles one finds thats 

I CRj . Ri.!) = U . I (Qi - Qi-1 ) ♦ v f ( Tl - T U1 ) ♦ S f (Pi - Pi.i) 
i=l i«l i=l i=i 

n 
■ - S £ (Wj, - W. ) 

i=l Ul 

K n -.*„ =U o C ^ -Q > ■♦ V T n - V* 
equation 18 for W n s 
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tf U oQn + V n ^ Mo ^ Vo p +w *n 

W n = s + P n - s ~ P * w o - T 

Substituting equation 15 in equation 19 one has J 

w %% + Vn ■ «n 

Equation 20 is accurate as it stands, i.e., no approximations have been made. 
It is equivalent to equation 5 only if j-g^j is negligible. Since equation 16 
tends to reduce |R n |, |R n | should eventually be sufficiently small, 'i.e., when 
the step has settled. When the step has settled, W is very nearly the desired 

sum. 

uv 
(2) PRODUCT. - The product W = ~- * P can be formed approximately by 

setting: 

AOj = AV t C21) 



ATi = AUj (22) 

W o = ~p* P o C23) 



AWj = sgn Rj sgn S (24) 

Again the sense of equation 24 is to reduce JR n j . Substituting into equation . 1 
we haves 

R i = K i-1 + u i Av i t Vi-iAUj + SAP A - SAWj 

= Ri.i ♦ "iAVi + v ul AUi + v^jUi.j . Uj^y^j ♦ sapj - sawj 

= R i-1 ♦ °i v i- u i-l v i-l ♦ SA Pi - SAWj 
= Ri.i «■ A CUV)i + SAPi - SAW t 
Summing - 

R n - «o = "nVn * SP n - U V . SP ♦ SW - SW n 

R n =Vn *SP n - SW n 

w n - s P « ' IT " ' *25) 
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Equation 25 is exact; no approximations have been made in its derivation. If 

R 
this step is settled, g"- should be small enough that it can be neglected and 

W n is approximately the desired product. 

(3) QUOTIENT. - As mentioned above, division is effected implicitly from 
multiplication. To find the quotient U = S ^\" P \ the product UV is compared 
to S(W - P) to determine whether U is too large or too small. This quotient 
can be formed by setting: 

AQi = AVi (26) 

ATi = AUi (27) 

Uo v o = S(W - P ) (28) 

AU i+1 = -(sgn V Q ) (sgn R.) (29) 

Substituting into equation 1 

R i = R i-1 + u i A v i + v i-l Au i + sA p i - SA w i (80) 

which is identical to the corresponding equation for multiplication. As for 

multiplication, summing equation 30 yields s 

Rn - R = u n v n ♦ SP„ - U V - SP Q + SW - SW„ 

Rn = U n V n + SP n - SW„ 

W n — P n R n 

u n = s V V s (31 > 

*n v n 

No approximations have been made in the derivation of equation 31, Therefore, 

R n 
U n is different from the desired quotient by the round off error —*, 

v n 

(4) SQUARE ROOT. - The square root algorism is effected by subtracting 
the product of the dependent variable times itself from the independent variable 
input. The sign of this difference indicates the next change in the dependent 
variable. The square root U ~«v/S(W - P) is formed by setting 

Uo = V (32) 

ATj = AQ 4 = AUj = AV A (33) 
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in 



U 2 = SW Q - SP (34) 

AU i+1 = -sgn Rj sgn V Q (35) 

Substituting in equation 1, 

R i = R i+1 + u i AU i 4 u i-a AU i + SAP i ' SAW i 
= R i4l + A(U 2 ) A + SAPj - SAW A 

Summing 

R n = U n 2 ♦■ SP„ - SW n - U 2 - SP Q + SW Q , 

U„ 2 = SW„ - SP„ ♦ R n (36) 

When R n settles to a negligible value, U„2 is approximately S(W n - P n ) , There- 
fore, U n must be close to the desired square root. Notice again the implicit 
approach to the extraction of the square rpot. 

(5) INTEGRALS. - Integrals can be approximated in the incremental com- 
puter in a single step in several ways. One way is to form sums of the inte- 
grands , That is, let U, V, P, Q, and T be independent variables. Solve equa- 
tion 1 for W. Then: 

R i = R i-1 + u i A Qi + v i-l AT i + SAP i - SAW i <*> 

2 Ri - Ri-1 = t (UjAQi + Vj jAtJ + S I (APi - AW.) 
i=l i=4 V i=l 

R n = £ UjAQi ♦ E Vi^ATi + SP n - SW n , SP * SW 

1=1 i=l 

Neglecting R„ 

IS ' 



1=1 i=l 

so Equation 37 is the Euler approximation to: 



(37) 



T, 



W n = s J UdQ ♦ | J VdT ♦ P n - P ♦ W 

On T„ 



(38) 
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A second way to integrate is the following* 
Let AUi = AV| 

AQi = AT 4 

Uo = V 

AW i+1 = sgn Rj[ sgn S 
Then equation 1 becomes 

R i = R i~l ♦ u i A <h ♦ u i-l A ^i * SAP i " SAW i 

Z Hi - Ri 1 = Z <«i ♦ U 4 ,) AQi + S I (AP A - AWj) 
i=l i=l " i=l 

R„ = Z (Oi ♦ U i . 1 )AQ 1 * S(P n . P - W n ♦ W ) (39) 

i=l 

Neglecting R 

«„ = S I (U i + "i-P A «i + P n - P o +W o M0 > 

Equation 40 is the trapezoidal approximation to: 

W„=| / "MQ + P n -P 0+ W o Ml) 

^0 

In those cases where the integrand U is an incrementally single valued function 
of the independent variable Q, the trapezoidal form is drift free, i.e., if the 
variable Q changes from its initial value and returns to its initial value, W n 
will become W Q if U is a single valued function of Q. An example of a single 
valued function is a polynomial in Q. An example of a function that is not an 
incrementally single valued function is sin Q which is obtained by double inte- 
gration rather than polynomial approximation. 

(6) DERIVATIVES. * The incremental computer is capable of taking time 
derivatives in a single step if properly stabilized. The restrictions of the 
algorism are: 
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AVj = APi = 

AQi =- ♦ 1 

A^i * ATj a - sgn Rj 
Substituting in equation 1: 

Ri * R ia ♦ Ui ♦ V AUi - SAWj 
Summing 

i=l 

Neglecting R n 

Z Ui ♦ V (U„ -*W =S(W n - W ) 



The term S Dj is approximately: 

1=1 .-. 

1 
18 ,' / t n + 2' 

2 Uj itf J . c U(T)dT 

1=1 ^ 

Then approximately: 

J n2 mr)dT + V (U n - U ) = sor-n - .Wo). 
2 



But OCt n «• 1) = U(t n ) ♦ ^UCt n ) + ^U(t n ) «■ |g"Ct ffl ) - -- - 

dW(t„) i dU(tJ i * i » ■ ■ 

So 0(t B ) -S^JL. (V p p^).^^.Jo..^O... 

■ ■ // /// 
sf = 0(t n ) ♦(»,♦&* *.||.*t... 



1,2 



Now UCt ra •♦ V ♦ i) = U(t n ) * CV ♦ ^) 0(t„) ♦ ° 2 o - -. 



So S |H = U(t„ + V + |) + 



I (V ° + £ 
8 2 



fi... 



S $ = "tt n ♦ V +%> - £ ( v + 1) u 



In other words U(t) is approximately equal to S times the derivative of W 
at time (t "♦ V + -i) . The V which was necessary for stability introduces a 
time lag. To make this lag small, one must try to make V small. The size of 
V is determined by the nature of W e If V is too small and the second deriva- 
tive of W is too high, U(t) will vary to both sides of the correct derivative 
value, eventually settling to the correct one. If V Q is too large, the delay 
will be unnecessarily long. 

(7) INTEGRAL WITH A RECIPROCAL INTEGRAND. - A variation of the integral 
algorism is the form: 

S /• dW 






(42) 



i 

2 



which is affected by letting; 
AT, = 



AV| = auj 
AP 4 = 
"o=V 



+1 = - sgn Rj sgn U Q 
Substituting in equation 1 

R i * R i-l ♦UiA^i +U i-l A Qi -SAW; 

R i - R i 1 + SAW, 

An - -J Ll! 1 

Summing over in cycles; 



C44) 
(45) 
(46) 
(47) 
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n o - s V Aw i + y R i - R i-l (48) 

If one neglects the contributions of the term containing the it's, equation 48 is 
approximately the integral equation 42. 

Through the use of integration with a reciprocal integrand, the logarithm 
can be generated. If we let Uj = Wj the form becomes; 

Q = 2) W *49) 

or = 2 log e W (50) 

(8) EXPONENTIAL. - The exponential is found from the logarithm algorisms 
by letting Q be the independent variable and taking the increments of W from the 
step as the dependent variable. The algorism is identical to that of the loga- 
rithm except that instead of equation 47 the dependent variable is determined 
bys 

AW i+1 = + (sgn R)(sgn S) 

e. SPECIAL ALGORISMS. - The following algorisms are special in the sense 
that they are not the usual arithmetic operations. They have been developed 
for some special purpose during the evolution or study of the fire control ap- 
plication. They are included here for completeness, as examples of the varia- 
tion of the incremental computer algorism, and for possible application in ser- 
vicing or other programs. This section can be skipped by the reader without 
losing the continuity of this presentation. 

1^ .Memory. Some time during each major cycle each position used in the random 
access memory has its old increment value removed and the updated increment 
entered. On some occasions it is desirable to have the previous value avail- 
able. It is possible to do this by using a step in the following way: 
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Set AUj = AV A =0 

4 = Vi =0 

5 = 1 

AW l+1 = sgn Rj 
That this is truly memory is proved as follows; With these restrictions equa- 
tion 1 becomes: 

«i =Ri-l * AP i - AWj 
If for i = n, R n = AP n - 1, then 

Vl = R n * ^ p n + l " Aw n 

= AP n .- 1 + AP n+1 - sgn ( P n - 1) 

= £P n - 1 ♦ AP n+1 - AP n 

= **„♦! -1 " 
since sgn (1 - 1) = sgn (0) = +• 1 
and sgn (-1 - 1) = sgn (-2) = - 1. 
but initially 

R = 0, and sgn R = + 1 = AW 

Rj = + AP 2 - 1 = APj - 1 

Therefore, equation 50 is proven by mathematical induction. In this case sgn 

CAP n+1 - 1) = AP n+1 = AW n<fr2 

or AW n = AP n .! 

In other words, on the nth major cycle, if A P n is desired, A P n is addressed 

directly; if AP„ , is desired, AW„ is addressed, 
" n-i ' n 

CI) SIQN. - Set 

"o = v o. A "i = A V is AT A = -AQ A 
APj = 0, AWj = sgn Rjj, S = 1 
Then Rj = R^j «• UjA Q A - U^jAQi - A W| 
= Ri-1 ♦ **l AQ 4 - AW 4 
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If for i = n 

R n =AU n AQ n -l, .(52) 

then AW n+1 = sgn (AU n AQ n - 1) = AU n AQ n (53) 

and R n+1 = AU n AQ n - 1 + AU n+1 AQ n+1 - AU n AQ n 

"n*l = Au n+1 A Qn+l ' i 
and since R =s 0, Rj = + AU^Qj - AWj = AUjAQj - 1, equation 52 is 

proved by induction. Therefore, equation 53, giving the desired product, is 

proved. 

(2) ABSOLUTE VALUE. - When a variable changes its magnitude by one, its 
absolute value changes by one. If the variable is positive, both changes are 
of the same sign; if negative, opposite sign. The change in the absolute value 
of a variable upon a change in the variable by one is given bys 

A ± |x( = sgn Xj sgn AXj (54) 

or A i|*i| = sgn X A AXj (55) 

One can form the product XjAXj by letting Uj = Xj and AQj = AXj, In order 
to sense the sign of this product, it should be put in a remainder alone. In 
other words a step which adds to the remainder the product X^Ax. and subtracts 
the previous product X. ,AX^ , has as its remainder X.Ax^, The signum func- 
tion of this remainder is the increment of the absolute value. 

We have seen above that an increment can be stored for an extra major cycle 
through the use of memory. One might expect a memory step would be required to 
have available the AX^ j required to form X._ jAX^ ^» * n one special case this 
memory step can be avoided. This case is the following; 

If a variable, X, is calculated in step j, the sign of its remainder is 
Stored in the random access memory on the fifth digit period after the start of 
step j ♦ 1. At any time prior to the fifth digit period the old value AX jj is 
available. After the fifth digit period the new value AXj is available. 
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During the step j ♦ 1 the increments for the calculation of step j ♦ 2 are be- 
ing selected. Step j + 2 can be used to calculate the absolute value of X| 
without a memory step by letting; 

AU X = AV A = AQj = AXj (56) 

ATj = -AX^j (57) 

U = V = ♦ 1 (58) 

S = (59) 

The accumulation of the sign of the remainder of step j + 2 is the absolute 

value of X e This step does not require servoing. 

To prove that sgn Rj is Aj X , consider the following: 
If equations 56, 57, and 58 are substituted in equation 1, it becomes; 
R t =R i . 1 +XJAXJ -Xi.jAXi.! 
Summing 

**n - **o = ^n^ x n - X AX 
R n = X n AX n - 1 
If |x n | > 2, sgn R n = sgn X„A X fl 

and sgn X n AX n = A n X since if X changes by one with the same sign as X it- 
self, its ibsolute value increases by one. If AX and X are of opposite sign, 
the absolute value decreases. 

If X n = + 1, sgn R n = AX n , since sgn = + 1 8 If X fl = ♦ 1, AX„ is A n |x|, 

If X n = 0, IajJ = 0, and ^ n |x| » - 1 since |-l| = Jl| =t + 1. In this case 
sgn R n = sgn - 1 = - 1 = A n |x|. 

Lastly, if X n = - 1 and A n X = + 1, sgn R n = sgn (-1-1) = - 1 = A n |x|. If 
X n =r - 1 and A X n = - 1, sgn R n = sgn (+ 1 - 1) = sgn = + 1 = A n |xL 

Therefore, for every value of X n , sgn R n = ^ n |x|. The accumulation of 
these increments is the absolute value of X„ 
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NOTEs The restriction that the calculation of the absolute value should 
occur two minor cycles subsequent to the calculation of the variable does not 
detract from its usefulness sine® this' timing is usually the most desired, 

(3) INDEPENDENT VARIABLE (TIME) CONIROIi, - In some servicing problems 
it is desirable to have an independent variable change from one value to another 
and remain within one Increment of 'that value Indefinitely, That is, the change 
in the variable will consist of N + i's followed by an indefinite sequence of 
alternately - 1 and ■'♦ 1, In order to perform this function three core memory 
positions are used. 

One memory position is used as a constant * 1, Any position not otherwise 
used can so serve, A second memory position Should be a "step function", i,e,, 
* 1 for the first major cycle and - 1 indefinitely thereafter. This latter 
function can be obtained by letting 

d -v »i-:-. 

AUj a AVj = APj = =s S 
AQi « . 1. 

AT U1 a", sgn Rj 
Then Ej a - 1 , ■- 1 a - 2 

' ■ Bn - B n __j - 1 ♦ 1 » R ra-i a .2 

That is after the first major cycle, sgn R m » - 1 - indefinitely. 

Given these two functions, the independent variable can be controlled in a . 

i ■•■■.■■■ 

single step by letting? 

■V s v o = N ■ " r 

AUj = AVj = AP| a o 
AQi = ♦ 1 
AWj a sgn Rj_j 



2e 
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S = 2 
and ATj = + sgn remainder of the step function. 
Then R x = N(+ 1) ♦ N(+ 1) - 2 = + 2N - 2 

R 2 = ♦ 2N - 2 ♦ N(+ 1) ♦ N(- 1) - 2 « ♦ 2N - 4 
and if n £ N + 1 

R n = R n _ 2 ♦ N - N - 2 =■ R n-2 - 2 =? ♦ 2N - 2n 

R N = 

R N+1 = - 2 

R N+2 = .2+N-N*2 = 

R N+2i = -2+N-N + 2=*0 

R N-2i+l = 0*N - N - 2 = - 2 
Therefore, the sign of this remainder gives the desired function. 

(4) GATING. - In some service problems it is desirable to gate certain 
variables. By "gating" is meant either allowing a sequence of increments to be 
duplicated exactly in a certain memory position until a certain time after which 
that position will have alternately plus ones and minus ones unconditionally, 
or causing a memory position to have alternate plus and minus ones until a 
certain time after which a given sequence of increments is duplicated. These 
gates will be referred to as "initially open" or "initially closed" respectively. 

Either of these gates requires a gate timing pulse. Incrementally the gate 
timing step memory position has alternately plus one and minus one until the 
switching time. At the switching time there will be one extra plus one and the 
alternate sequence will continue thereafter. To effect such a function the 



*• algorism is: 



«0 = 



V =3 
S = 2M 
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AUi = -AV A = AWj = sgn R^j 
APj = 
ATj = -AQj 
We wish to have this timing sequence switch when T N = M, Let T be a .'..variable 

such that T = and T>0, To demonstrate that this gate timing step produces 
the desired sequence, substitute in equation 1, Then 

■Rj sQU C-ATj) ♦ 3ATi - 2M = 2ATj - 2M<0 

R 2 - 2AT i - 2M ' oAT 2 * 2AT 2 * 2M ~ 2 ^ AT i * At 2> = 2 CT2 - Tq) £ 
As long as the sequence is alternate 0| will be ♦ 1 for odd major cycles and 
for even major cycles. Likewise, V*n, will be * 3 for odd and '♦ 2 for ©vera 
major cycles^ If we assume that R 2 i^® where i is an integer, then 

R 2i*I * R 2i - AT 2i+l * 3AT 2i*I * 2n = R 2i ~ 2M * 2AT 2i*l 
and if R 2i ~ 2M ♦ 2^T 2i *2< 0. 

R 2i *2 = R 2i- 2M * 2AT 2iUl ' ° * 2 AT 2i*2 = R 2i * 2CT 2i*2 - T 2i> 
If we assume the sequence of alternate increments, the general values of the 

remainders are; 

»2n = 2W 2n - V =2T 2n 

'«2«*1 3 »2n*l " 2M \ 

Conversely, the sequence of sgn R* will be alternate as long as sgn R 2 n*l 1* 

- 1 where i£2n ♦ 1. The smallest value of i for which this condition does not 
hold is found as foliowss 
' 2T 2n+1 - 2M^0 

T 2n*l^ 
In other words, the alternate sequence will continue until the first T on an 

odd major cycle is greater or equal to flfl, the predetermined switching point. 

Assume that switching occurs on the Nth major cycle, i„e eo T N ^M 

Then R N = 2T N ^ * 2AT N - 2M = or 1 



22 



But sgn = sgn 1 = + 1„ The sequence sgn Rj is no longer alternate but con- 
tains a + 1 for the Nth major cycle as well as for the N-lth, Then U^ a 2, 
The next remainder is; 

R N+1 = R N . 2AT N+1 * 2AT N+1 - 2M = - 2M or 1 - 2M 

V N+i = io Since i ~ 2M<0 » U N+2 = lo Usin ° the8e values 

R N+2 = R N - 2M - A T N+2 ♦ AT N+2 * 211 = R N » or 1 

After r the N ♦ 2th major cycle U N+2 = 1, V N+2 = 2 * and R N+2 = 0. These are 
precisely the values obtained after the Nth major cycle. In other words, the 
step has become cyclic with a period of two major cycles . The remainder of 
this step is alternately and - 2M Cor 1 and 1 - 2M), In either case sgn % 
will be an infinitely long alternating sequence. Therefore, this step yields 
the desired gate timing pulse. This step is not reversible; the sequence will 
remain alternating regardless of subsequent values of T, 

A single gate timing step can supply timing to open and/or close any number 
of gates at some particular time. Both the initially open and the initially 
closed gates will be considered here. 

The initially opened gate is formed by-^tke •foilwing algorisms 

Vo = 2 

S = 2 

AQj = - AT t 

AUj a -AVj = sgn Rj in the gate timing step, 

AP 4 = 

AWj = sgn R| i 
Prior to switching the Uj will be 1 on even major cycles and on odd, Vj i 
will be 3 on even major cycles and 2 on odd. In every major cycle prior to 
switching V^j - U i = 2 C Substituting in equation 1 we haves 
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Ri = R. x + UjC- ATj) .♦ V A j( ATj) - 2AW i 
= Rj.i + 2ATj - 2AW A 

R^ = ♦ 2ATj - 2 
But sgn Rj = A Tj 

So R 2 = 2AT X - 2 «• 2AT 2 - 2ATj = 2AT 2 - 2 

In general 

R n = 2AT n - 2 
and sgn R n = A T n 

prior to the gate switching. If the gate switching pulse occurs on the Nth 
major cycle where N is odd, the values of U and V will have the following 
values in the neighborhood of the Nth major cycle. 

Major Cycle si Uj Vj V A-1 v i-l - u i 
N - 3 1 2 3 2 

N-2032 2 

N - 1 1 / 2 - 3 2 

. N 2 12 

N '♦ 1 1 2 1 

N ♦ 2 2 12 

After switchings 

R N = 2AT N j - 2 - 2AT N + 2AT N - 2AT N _j = - 2 ■' 
n ml = -2*2 = 
On ail subsequent major cycles ATj has no effect because Vj_j - Oj = and 

k 

T the sgn Rj will alternate^ Therefore, if k£ N, £ sgn Rj »Tjj i - T k even 

0/5 i=l . 

8 = Vl - T - 1 k odd, 

The initially closed gate is formed by letting the algorism bes 
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Un = 



V = i 



S = 2 

AQi = - AT| 

AV| = -AUj = 8 9 n R i * n the gate timing step 

AP ~ 

AWj = sgn Rj_| 
Prior to switching both the Uj and the V|_| will be 1 on even major cycles 
and on the odd. The difference Uj - V^j = ©„ Substituting in equation Is 

*i = *i-l * °i A Qi - V i-1 AQ i - 2AW £ 
= R i-pl ♦ OAQ* - 2AWj x R^j . 2AWj 

Rj s o - 2 = ■- 2 

R2 = - 2 *.2 = 

R 2n = 

ft 2n*l = " 2 

If the gate switching pulse occurs on the Nth major cycle where N is odd, the 

values of and V will be the following in the neighborhood of the Nth major 
cycles 



Major cycle » i 


«v 


v i 


¥ i-l 


°i - ¥ i-i 


N - 3 


1 





1 





N - 2 





1 








N - 1 


i 


o • 


1 


■0 


N 


2 


-1 





2 


N ♦ 1. 


1 





-1 


2 


N*2 


2 


-1 





2 


Rflj = * 2AQ N - 


2 








Sgn R N = A Q N 











2S 



R N*1 = 2 AQ N - 2 ♦ 2AQ^ - 2AQ N - 2AQ N+1 - 2 
sgn R N+1 a AQ N+1 



In general 



s«n r n+j = q n+j 

The function 

k 

X sgn R A = Q k - Q N where k £ N 

which is the desired gated function. 

(5) ERRORS IN INCREMENTAL ALGORISMS . - In performing the incremental 
algorism the computer can come up with only one possible answer. Since all 
the operations are digital, there is not the slightest ambiguity about the 
answer. This answer will be the same whether the incremental algorism is per- 
formed by a small, specially built computer or a large general purpose computer, 
However, the function generated by the incremental computer may not coincide 
with the function for which its program was designed. The difference between 
the desired answer and the answer the incremental computer will yield with a 
given program is called the program error. 

In any digital computer there exists round off error when a number is 
handled whose significance exceeds the ability of the machine to represent the 
number as an integer. For example, a ten binary digit number can represent any 
integer from one to 1024. If ten binary digits are used to represent numbers 
fromO.Ol to 10.24, the number 7.8382 would have to be rounded off to 7.84 with 
a round off error of 0.0018 or 0.01. 
"f In an incremental computer the round off error is usually no more than plus 

!/> 

K or minus one increment on input quantities and slightly more on output quanti- 

se 

ties, assuming that the computer can keep up with their changes. In order to 

keep up with the changes in the variables, the increments must be at least as 
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large as the maximum changes of the variables that can occur in one major cycle. 
For example, the error of a variable input or output due to round off is about 
equal to the maximum change the variable can make in 0.005 seconds in the case 
of real time application. In many applications this amount of error is small 

Round off error in an incremental computer corresponds to the error made in 
an analog computer due to the difference between the analog representation of 
a number and the number itself. To decrease the round off error by a factor 
of two in an analog machine may be a very difficult engineering operation. To 
decrease the round off error by a factor of two in an incremental computer not 
used for real time control, one more binary digit is required for the numbers 
and twice as much time is required for the computation. 

The round-off error for the outputs from addition, subtraction, muitiplica- 

R 
tion, and integration steps is equal to the term g. This term is usually less 

than one, but it can be as high as three even when the numerical process has 
"settled". In order to keep the fractional error small in addition, subtrac- 
tion, multiplication and integration, the answer should be kept large as com- 
pared to the possible error. That is, the magnitude of the quantities 

U Q + V T nv 

-2 — - — °— + P for addition and subtraction, ~ ♦ P for multiplication, and 

^jUdT for integration must be kept as large as is consistent with the require- 
ment that the change in the quantity shall not exceed one in one major cycle. 
The only quantity which may be adjusted to accomplish this is the scale factor. 
Therefore, to minimize round-off error in addition, subtraction, multiplication 
and integration, S should be just large enough that the arithmetic process can 
T keep up with the changes that can occur. In division and differentiation the 

x quantities to be kept large ares * W „" ^ and S^, In these cases S should be 

gu v at 

made as large as the requirement that the step must keep up will allow. Round- 
off is the only error that occurs in addition, subtraction, multiplication and 
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division. The round-off term is not discarded as in normal round-off, but is 
retained, preventing the error from accumulating. The error in the answer is 
due only to the final rounding off. If the independent variables are returned 

to their Initial values, the dependent variable will return to its initial 
value. 

Round-off error is not the only error that occurs in incremental integra- 
tion and arithmetic processes which use integration,, More accurately, integra- 
tion accumulates the round-off error which may grow systematically or statisti- 
cally. The integrated round-off error is called drift. To investigate drift, 
consider the polygonal function: 



Then 



FCT) = FCTi) ♦ CT - Ti ) fr ( V ^*+ l > - FCTi) 1 



C60) 



for Tj£T£Tj +1 
Ti>T>T i+1 

and T i+1 = Ti ♦ AT i+1 



if AT i+1 = * 1 
if AT M = - 1 



/%CT)dT= / Ti W ) ♦ [T - Tj [^^ 
*2 



l l i 2 l J T n-1 



F<Vl ,„„ . x„,, \ '"l ; ZV' } 



- TT 



n-1 



•n-1 



28 






T i * T ' 

-*- - T,T ♦ — 8- 
2 A l i o 2 



2 



= FCT )(Tj - T ) ♦ [FCrp - F(T )] 

[ T n* TT la-lf] 



2 



)CT,-Tj*ct,.t > (g » ) : r M .. 



'ii *" x n-l 



*F(T )(Ti - T ] 



1 " "0 



T. - T. 



* CT n - T n .!> F(T„.!) ♦ " ~ 2 n -t [F(T n ) - F(T„.i)] 



-CI.-V £Hl>4£<V ♦ ... ♦cr.-vp 



■ FCTi) + F(Ti.i) 



F(T n ) + r(T B ,i) ' 
2 



That iss 



i=l 



T. 



2 



F(T)dT 



i=i 2 * 



is suggestive of the incremental integration formula, equation 39, 
equation 61 into equation 39, 



- "o - 5 



UdT - -f 
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If U is the polygonal function: 



+ (T - t|) 



UCT i+1 ) - U(Ti>" 



AT i+1 
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T<< T<T 



i+1 



is + 



R 

If the process has settled in n major cycles, *p can be neglected and 

T 
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Equation 64 is accurate except for the round-off error ~* for the polygonal 
function defined by equation 63, In many applications, however, the function 
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to be integrated is not of this form. Let the function to be integrated be 
G(T). Let the changes of G from one major cycle to the next have a magnitude 
no greater than one. The polygonal function U(T) can be made to approximate 
G(T) within one increment. Define this difference ass 

S(T) = U(T) - G(T) (65) 

S(T) is the input round-off error. 

Then; I J " G(T)dT = | / " 0(T)dT + I / " S(T)dT (66) 

T S T T 

2 rn 
The term o J S(T)dT represents the error that would be caused by the approxi 

T ■ 

I T 

S(T)j * 1, J S(T) 

I T 



dT 



£ n. Moreover, the 



actual error is usually much less than n because the average value of S(T) is 
approximately one-half in most cases instead of one, and the sign of S(T) tends 
to alternate. If one assumes that this error accumulated statistically and 
that the probable error for the round off error for one major cycle is or f then 
the probable error for n major cycles is WrT. Since the probable error a 

2 f T i+l 
is the probable error of the term § J S(T)dT, its value for most appli- 

cations is about 0.0005, Under these circumstances an accumulated probable 
error will be equal to one increment in 4,000,000 major cycles (5.6 hours of 
computer operation). 

To show that T is not a single valued function of W, it is sufficient to 
show that for two different major cycles with the same value of W there are 
different values for T. Consider a sequence of AW's such that they are ♦ 1 
until W n is reached, and are - 1 immediately thereafter. In this case .... 
AW n = 1, Aw n+1 » - 1. Then W n+1 = \H q1 + Aw n + AW fl+1 = W^. Without loss 
of generality assume that 
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0<R n-2 <s * 2. (67) 

Assuming that W is positive for this problem, equation 57 

AT i+l = - sgn W sgn Rj (57) 

indicates that A/T , = - 1 # Substituting in equation 41 

R n-1 - R n-2 + < w n-l + W n-2> AT n-l " SA Vl 
= R n-2 - 2w n-2 - 1 - S <0 
Therefore, 

AT n = + 1. 
Then R n = R n _ x + (W n + W nl ) AT n - S AW n 

= R n-2 - 2W n _ 2 - 1 - S + 2W n _ 2 + 3 - S 

R n = R n _ 2 + 2 - 2S<0 
Therefore, 

T n+1 = + 1. 
Then T n+1 5 T B-1 + AT n ♦ AT n+1 = T n-1 + 2 £ T nA 

Although W n+ j = W n _^, T n+ j ^ T n _j # T is not a single valued function of W. 
Therefore, the drift in the logarithm is not reversible,, 

(6) SUMMARY. - The algorisms of addition, subtraction, multiplication, 
and division are subject only to round-off error. These errors are small, e.g., 
with full scale numbers equal to 1000 the round-off error is about 0.1 percent. 
The algorisms involving integration are subject to drift as well as round-off 
errors. Drift is either reversible or irreversible and it is either systematic 
or random. (Actually, there are no random processes in an incremental computer. 
By random is meant that which cannot be predicted analytically from the dif- 
ference equations describing the process.) The random drift is very small and 
can usually be neglected. The systematic drift can be anticipated analytically, 
and perhaps eliminated by a suitable modification of the algorism. Reversible 
drift occurs when a single valued function is integrated. Such drift is zero 
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when the independent variable is returned to its initial value. The algorism 
for differentiation is unstable, but can be used with suitable feedback. 

By using the smallest value of a scale factor that allows the step to keep 
up in the algorisms for addition, subtraction, multiplication, and integration 
and the largest scale factor;for division and differentiation, the round-off 
errors can be kept to a minimum. 

f e APPLICATIONS OF THE INCREMENTAL TECHNIQUES USING SEVERAL STEPS. - In 
the following discussion the steps for the variables are indicated by a super- 
script. As before, the number of the major cycle is indicated by a subscript. 
Independent variables will be indicated by letters without superscripts. 

(1) FILTERS. - A single RC low pass filter satisfies the following 
differential equation: 

fe j - e o = BC JT 

T 
or RCe = / (ej - e )dt (68) 

By using a single integration step with two simultaneous point slope (or Euler) 
integrations, a single stage RC filter represented by equation 68 can be 
realized. The algorism for this step is: 

U = ( ej ) AUj = (Ae^ 

AQ A = - ATj = + 1 APj =0 

V = AV A = AWj = sgn R U1 

S =r RC ♦ 2 
This algorism corresponds to the integral equations 

8 r r 

x J e.dr -J e d = (RC ♦ 2)e 

f f 

or J ejd -J e d » RCe (69) 
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which is the filter differential equation 68. Therefore, when equation 69 is 
satisfied Wj* .is very nearly e „ These filters can be cascaded to give rather 
sharp attenuation. By cascading five sections, the mathematics simulates five 
electrical RC filters cascaded with unity gain isolation amplifiers between sec- 
tions* The circuit would look like; 




s program was run numerically on a simulation routine in an 1103 large scale 
computer to obtain the frequency response of the system and the five outputs 
agreed in amplitude and phase shift with that expected from filter theory, 

(2) POLYNOMIAL. - Since multiplication can be performed in a single 

o 
step, a quadratic can befformed in one step. If it is desired to form Ax** ♦ BX' 

* C, form the product 

W n - W = AxCx «• |) and let W =? G, 

Then W is the desired quadratic. In specific incremental notation let 

D o = V v o = x o + v AP i = °. AU i = AV i = AT i = 

B 
signum Rj, Then U = x, V = x + j, and 



I = Ax i9 an« 



a*j 



x(x + p 



when the step has settled down. 



The scale factor, S, is chosen as a compromise as indicated by the following 
discussion. For the proper choice of the initial value W ,' W is proportional 

to the desired quadratic. A cubic can be generated in two steps. Probably 

3 2 
simplest way of generating the cubic Ax + Bx + Cx + D is to generate the 
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quadratic Ax + Bx + C in one step, and multiply it by x and add D in the sec- 
ond step. A more general way is to factor the cubic into 

(x ■*• a) {ax 2 + (B - aA)x + (C + a 2 A - Ba)} ♦ D - a 3 A ♦ a 2 B - ac (70) 
The quadratic indicated is generated in the first step and the quadratic is 
multiplied by x + a and the constant term is added to it in the second step. 
Factoring the quadratic in this manner is helpful because the usual limitation 
on the accuracy of such a computation is due to the scaling of the quadratic 
which is necessary in the first step to allow the step to keep up with changes 
that occur in the independent variable. This rate of change is given bys 

i-JAx 2 + (B - aA)x + (C + a 2 A - Ba)) = (2Ax + B - aA)|4rl C71) 

dt l J ,dtr 

! dxi 
where hrr is the maximum absolute value of the rate of change of x at this 

value of x. The constant "a" is chosen such that the maximum rate at which the 

quadratic can change at all values of x is as small as possible. If l-j*-! is 

constant throughout the range of x v equation 71 is linear in x and a = j + x,^ 

+ x max . For any choice of a, if the scale factor S is equal to the maximum 

rate of change of the quadratic, the quadratic will keep up with the changes 

of x throughout its range. For this choice of a, the quadratic will keep up 

with the minimum of round-off error. In the case of a constant |~|. 

Idtl 

S = A * x max - x min* Idtl' 

Quartics and higher degree polynomials can be generated as well. In gene- 
ral a polynomial of a degree n can be generated in n-1 steps. Variation in the 
details of generation can be made allowing the programmer some flexibility in 
adapting the computer to a particular polynomial. For example, a sixth degree 
polynomial can be formed by the product of two cubics, the product of a quartic 
and a quadratic, the product of three quadratics, or the five steps can produce 
the sequence of quadratic, cubic, quartic, fifth degree, and sixth degree. One 
of these systems would be the best for a given sixth degree polynomial. The 
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an 



a. 



best system would be the one which can keep up with the change in the indepen- 
dent variable with the minimum of error in the polynomial caused by round-off 
of all the steps. To obtain the various constants required for the optimum sys- 
tem, criteria of the type suggested for the cubic should be used. These cri- 
teria are too complicated to be included in a paper of this scope, 

.2 
(3) SINE AND COSINE. - The differential equation ^-4 + aj 2 y = (72) 

dx* 

defines the trigonometric sine and cosine. This equation can be solved by the 
incremental computer either through the use of the differentiation routine or 
through the integration routine. Inasmuch as differentiation requires that the 
independent variable increments always have the same sign, integration is the 
more useful. Integrating equation 72 twice we haves 

y Cx) = -^2 J dt { J d7}y(rj) ♦ cj «• C 2 C73) 

o 

The general solution to either equation 72 or equation 73 is: 

y = A cos ^ x + B sin u)x (74) 

If we wish to find the cos ui x, the initial conditions are y(o) = A. Then 

y = A cos wx, The first integration of A cos «>x yields A sin wx 5 The seejnd 

integration yields -A cos ujx. The sum of A cos a) x and -A cos m x should equal 

0„ The sign of the sum of the two calculated terms determines the sign of the 

correction to be applied to A cos '*»> x. 

The sine and cosine can be generated by the direct application of two steps 

» 

of integration. This algorism iss 
AP 1 (1) = A?/ 2 * = 



AT/ 1 * » A Qi (1) ..AT^ 2 ) =-AQ/ 2 > = A Xi 



x °o = V = A cos wx o ^77) 



U o C2> = ¥ o C2) =Asin^x ♦ 1 
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(78) 

AW.^ = sgn k/^ } (79) 



AUi (2) = AVi (2) = Aw{Vl = s 9 n R i U) (80 > 

S (1) = S (2) = |- > 2A (81) 

One incremental change in x represents an angular change of 4 radians. With 
this algorism the sine and cosine are generated very well under favorable cir- 
cumstances. This program was tried with several values of A and S to hcheck on 
its usefulness. For example, it was found that with A = 1000 and S = 2000, 
that after 6,283 cycles representing 2 it radians of angular change, the error 
in the function A cos 2 it was between two and three increments, about 0,25 per- 
cent. 

<, . ■ 

If, for some reason, it is known that the angle for which the sine and/or 
cosine is needed is limited by the problem, it is advantageous to change the 
scale factors so that they are no longer equal and are just small enough to 
allow both steps of the algorism to keep up with the change of the independent 
variable. If, for example, it is known that -20° £ -tu x £ ♦ 20°, we know that 

I dx cos <jUX l = I"* 8 * 11 tux I - w ° 342 
If ot = o.OOl, we can generate 2924 cos at x in the second step instead of 1000 

cos di x by letting the scale factor in the second step be 684 instead of 2000 
and the scale factor in the first step be 5848 instead of 2000, The variable 
1000 sin w x is computed in the first step just as before. By adjusting these 
scale factors the rounding off error is reduced for the cosine in this applica- 
tion causing less final error in both the sine and cosine. 

By looking more closely at the difference equation actually solved by the 
algorism, the source of this small error can be found. Notice that in equation 
79, A u/ 1) = AVj* 1 * = sgn R^ ^ In this case there is a lag, i.e., 

cos a> xj_2 + cos itt xj ] 
Aj sin ojx = ~ = - Axj 



sin (O Xi i + sin wxj 

whereas A{ cos cu x = — Axj. 

x S 

The latter represents trapezoidal integration. The difference equation actually 
being solved is: 



i-1 

Yi - y = - Z 
k=l 



k r i 

Zjyj * Xj-i] % + 




S 2 

b. 


* 



>AT k (82) 



which corresponds to the integral equations 

J dr)y(r)) +C X 



f X ~ l r 

(x) = -u> 2 J dr 

o 



♦ C. 



C83) 



or to the differential equation: 

i&£0 + ci 2 y(x-l) = 
dx 2 



C84) 



C85) 



Expanding equation 84 in a Taylor's series, 

s?d£y. + v(x) Msl + 1 d 2 Y (x) l d 3 Y <x) t 

4 dx 2 + y(x) " dx + 2 dx 2 "6 dx 3 ♦— ° 

This series is strongly convergent. An extremely small error is made by round- 

ing off equation 85 to: 

S 2 " 
4 



y + y-y + o =0 



C86) 



Equation 86 has the general solution: 



2x r 



[x) = e s2+2 A sin ^^S x ♦ B cos l£^l 
L 2 ♦ S 2 2 + S 2 



C87) 



In the usual case, S is of the order of 1000. With such an S, it is justified 

to approximate the solution by: 

— ■ <* 
y(x) = e [A sin g x' + B cos ^ xj C88) 

Equation 88 represents the solution to the difference equation 82 actually be*. 
ing solved by the algorism of equations 75 through equation 81. For many appli- 
cations this solution is close enough to the sine and cosine to be used. In 
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particular, in the usual case where the independent variable changes relatively 
slowly, the scale factor can be made large, and the exponential multiplier in 
equation 88 changes from one slowly. In the cases where S cannot be made large 
and the computation is to be carried on for a number of seconds, this algorism 
produces considerable error. 

For a numerical check an incremental algorism of this kind was computed on 
an 1103 large scale digital computer. All of the steps were performed to simu- . 
late the incremental operation. The initial values were A = 150, B ?0, and 
x = e In order to show how the algorism would differ from the .sine in a rea- 
sonable number of cycles, a small scale factor was required. The scale factor 
was 400, After 23,562 major cycles (corresponding to two minutes of incremental 

computer running time), y was found to be -200,9. From equation 88, 

2 i 23.562 
e 160,000 , 150 . (-0.9987) = -201. 1; in very good agreement. In this case, 

after a long time (two minutes), and under poor conditions (small S and no 

changes in AT) ( a 35 percent error was produced. Compare this result with the 

more realistic result above. 

In the derivation of difference equation 82 it was assumed that AT was ♦ 1. 

Since x is an independent variable, it is conceivable that AT can be both 

signs. During each major cycle that AT is different from the previous major 

cycle, the difference equation 82 does not apply, and there is no drift. 

Defining; 

N e no„ of major cycles. 

n s no. of changes in T. 
the equation 88 becomes s 
2(N-n) 



C" If" Q O "B 

y(x)=e I A sin g x + B cos c . x 1 



I 

o 



The variable n can be large when x is varying slowly causing the changes in x 
to be alternately plus and minus one. 

It is possible to alter the algorism to generate the trigonometric func- 
tions more accurately. 

In some applications x is monotonic, that is, AT = Ax = + 1 for every 
major cycle. Such would be the case, for example, if a sine function of a cer- 
tain frequency (relative to the speed of the drum) were desired. The sine- 
cosine algorism for monotonic independent variable can be formed by modifying 
the algorism to read: 

AT/ 1 * = AT/ 2) = APi (1) ="AP/2) = o 
AQi (1) = -AQi (2) = Ax. = + 1 
U = A cos aix 

U * 2) = A sin oj (x + j:) + 1 

s (l) =s (2) i ■ A 

u) — 

AD| (1) = A Wl <2) = S g„ B^J '■ 

AO/ 2 ' = Awft =sgnR i (l) , 

This modified algorism, using point slope instead of trapezoidal integration, 
corresponds to the incremental equation; 

D <»., . v 1 u J (2)Ax J = . y 1 ±H y -^1 u. (1) 

1=0 S< 2 > i= 0S (2) jtosd) J 

i=0 j=0 J 
In general y(t*l) - y(t-|) = y(t) + iy(t) + ~y + i-y + ^y + ..... 



- y(t) -^y(t) + |y(t) -^y + { ^-y + ..... 



= y(t) = & * ^ ♦ 



• • • • 
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y(t) = yU^) T y(t-J) - ^y - jfeg, * ..... 

I y(t) = y(t>4) - y(a^) - I (£jy ♦ jfep ♦ 
t=a t=a 



or approximately 



b r b+« f b+9 #/ f b^ v 

iy(t)= i 2 yHr)dr .1^ J 2*t--1- ' 2 yd r.... f 
t-a b-^ a-^ a- 2 

- .1 " ■ i " i *' i in i 

« / 2 v . {T)dr 7 (b ^ ~ ' (a "^ y <b ' ^ " y (b ^ , 

1 y ' " 24 1920 



a- 
2 

Then 



* (1) f "^ ' U'(i4) u'"(i4) 

Z Uj UJ = J U(r) d T _w -X .w—^L. ♦ .... + Const, 

j^O J a4 24 192 ° 



So 

n-1 i _ _ t n-1+1 



U) 



n-l i ,. v « f n-l+1 f n ' ., 41 

I S0 J CD S .«2/ d^/ -dTUCT) -«4 0.«i c 

i=0 j=0 ■ -1 -1 ^ 4 1920 



That is, 

CD 



r V dt /'dru^-^m-.... 

= **>* ./ 



-U 
n 



Differentiating twice 

(1 - m 4 ) U n <D ♦ a* 2 U n (1 > =0 

1 - a; * 

9 (j 2 

Let n * = -^ (90) 

1 - ut H 

Substituting in equation 90 

y Cx) . + Q, 2 y(x) = (91) 
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This differential equation is of the form of equation 72.; Equation 91 has the 
general solution y(x) = A sin £lx ♦ B cos fix. The approximations necessary 
to derive equation 91 from the integral equation describing the effect of the 
modified algorism are much less severe than those necessary for the unmodified 
algorism. The modified algorism is much more accurate than the unmodified,, . 
virtually eliminating the exponential term of equation 88 found for the unmodi- 
fied algorism* The only other errors of this algorism are the accumulation, of ' 
the integral of the difference between the sine and cosine functions and their 
polygonal representations and the final round-off error. It must be remembered 
however, that this modified algorism is most accurate when the independent 
variable is monotonic. With non-monotone variables, this method is more accu°» 
rate than the method using trapezoidal integration. 

To demonstrate these methods, simulations were carried out. In the case of 
a monotone variable, the two methods were used to generate 120 sin oj t. The 
simulations were carried out for 13,000 major cycles (65 seconds of real time 
and 13 complete 360° rotations of the input angle). Here again the scale factor 
was chosen low in order to show the errors in a bad case. The results were; 
the first method (trapezoidal integration) yielded 157.4 and the second method 
(point slope integration) yielded 120.38, Their relative errors were 31 perceat 
and 0.3 percent in this extreme case. 

A second simulation was done where the input variable was a sinusoidal fun®. 
'tion of time instead of monotone. The errors made by the second method were 
©bout one-third that of the first. 

• In any of the systems discussed, if w is very large, S must be made small 
in the" integration steps and the amplitude, A, will have to be small. This com- 
bination tends to yield more drift due to the differences between the desired 
trigonometric function and its polygonal representation. Both of these systems 
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are subject to this drift to the same extent, Since polynomials are not sub- 
ject to drift error, a polynomial approximation to the sine or cosine can be 
made which will give sufficient accuracy with final round off being the only 
computation error. If the independent variable is restricted to a suitable 
range, e.g., - tt < a; x £ ♦ 7r , the approximation can be made directly by. 
Legendre or Tchebycheff polynomials. If, j»n the other hand, the independent • 
variable has an extended range, the polynomial should be chosen to approximate 
the trigonometric function from ^ip- - - -V-i anc ' * e t 9 = ^x + 2mr 
where n is a positive or negative integer. The polynomial approximation to 
the sines 

sin x » .986x - . 143x 3 
is in error no more than 0.006 at any value of x between ~w- and : ~. This 
cubic approximation should be close enough for those cases where the use of 
the polynomial is called for. In the extreme case, where a very accurate sine 
is needed after a very large number of major cycles, it may be necessary to use 
a fifth degree approximation. 

(4) CONCLUSION. - To generate the sine and cosine, the simplest and 
generally most useful method is to use two steps for double integration. This 
method lends itself to a wide range of applications. In some cases this method 
may lead to a drift error that is excessive. 

A second method has been devised which greatly reduces this drift error in 
those cases where the independent variable is monotone. In every case of a 
monotonic independent variable this method should be used since it, too, re- 
<«* quires only two steps of integration, but yields greater accuracy than the 
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first method. 

In some cases, generation of the sine and cosine by double integration may 
not be satisfactory. Suitable approximations can be made with polynomials 
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either with one polynomial for the whole range or with a polynomial approxima- 
tion from z t£- to ^5- repeated over the entire range. 



t 
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2, PROGRAMMING THE INCREMENTAL COMPUTER 

a, INTRODUCTION. - This section is concerned with the preparation of a 
for computation on an incremental computer. This preparation may be 

.trarily divided into a series of interrelated operations. It is assumed 
that, as is usually the case, the problem is to generate some output variable, 
which is a function of the input variables. The first step is to determine a 
series of operations which will generate the desired function. This procedure 
is known as "programming". The selection of the proper operations and inter- 
mediate functions is subject to various criteria to be discussed below. Further, 
the optimum sequence of operations must be determined. 

The second step is to determine coefficients for each of the intermediate 
functions, a procedure known as "scaling". The scaling of a function determines 
the precision to which it is known in the computer and also the rate at which 
the computer quantity can vary. The accuracy of the output function is quite 
dependent upon this scaling operation. As the incremental computer computes a 
function relative to its initial value, the initial values of all the functions 
must be computed and the appropriate quantities inserted in each register. In 
cases where the inputs start out at values different from the initial values in 
the computer, it takes a period of time for the functions in the computer to 
equate themselves with the inputs. This period, known as "settling time" is 
kept to a minimum by proper choice of initial values. 

The various steps listed separately above are actually interrelated. In 
practice this means that the changes in one step affect the others so that it 
is often necessary to iterate the steps. 

Finally, the operations must be coded for the computer, and the input' tapes' 
must be cut. 



b, DEFINITIONS, - At this point some of the terms common to this work will 
pe def ined e 
Coefficient -This is a-' constant, which a function Is .multiplied -by, in the 

computer, A function x is represented by Cx increments in the computer, 

where C x is the coefficient of the function, 
A direct operation is defined as one in which AW is the dependent variable* 

An inverse operation is one in which A W is an independent variable. 
The output , or result of step, refers to either the sign of the remainder ©f 

the step' or else the dependent variable generated by the step and defined 

by the equation 

n 

*n - * * 2 sgn Hi 

1=1 

Overload ing- A step is said to be overloaded if the output function changes 
more rapidly than one increment per cycle. This condition is characterized 
by large values on- the R-line, It is the opposite of keeping uj e 

Scaling is the operation of determining scale factors (and/or coefficients) for 
the steps of a program,, 

Scaling down is to decrease the coefficient of the output of a step, scaling uip 
to increase it„ 

Settling - Immediately after a computation -has -started,- -many steps and inputs 
will be overloaded. When no more steps or inputs are overloaded,, they are 
ail settled. The settling U,me is the time from the start of a computation 
until ail steps and inputs are settled. 

Flipping a step is to reverse the sign of its dependent variable. 
c e PROGRAMMING. . - 

CD CHOICE OB* FUNCTIONS. - In preparing a problem for an Incremental 

computer the first operation is to write the program. This consists of a 
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sequence of intermediate functions which bridge the gap from the input functions 
to the output functions. In order to obtain the greatest possible accuracy in 
the output, it is necessary to choose terms of the sequence which introduce a 
minimum amount of error, and, also, to arrange them in an optimum order. 

In choosing a set of intermediate functions, it is desirable to select 
functions which fall in the general category of analytically "well-behaved" 
functions; that is to say that they should be continuous and reasonably bounded. 
Furthermore, they should be relatively insensitive to absolute error and should 
tend to minimize its propogation. For example, the relative error in the dif- 
ference between two almost equal terms is quite sensitive to variations in 
either one of them. Or, to cite another case, the square root of a function 
becomes highly sensitive to absolute error in the function as it approaches 
zero. 

In addition to the restrictions placed on the function from the general 
standpoint of error, the incremental computer itself has properties which 
should be taken into account. Due to the representation of a function which 
is used in the computer, absolute rather than relative error becomes the most 
applicable type. In the many cases where a certain relative error is aimed 
for it is usually desirable to select functions with as limited a range as 
possible so that the two types of error may be as similar as possible. This 
limiting of range serves other useful functions as well, for as the maximum 
rate of change of any computer variable is one increment per cycle, a limited 
range allows a much greater total number of increments to be used to represent 
the function, other things being equal, for the same settling time. Scaling of 
a restricted function is usually more satisfactory as wide ranges tend to allow 
scaling to be too high at one end and too low at the other. This restriction 
of functions applies not only to the functions themselves but to their rates as 
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well. Here, even more clearly, the accuracy with which a function may be repre- 
sented is inversely proportional to its maximum rate. This may be seen from i 
the fact that the round-off error and the maximum rate are both proportional. to 
increment size. 

The principles which make inverse pairs, e,g, , multiplication and division, 
s© easily obtainable in the Incremental computer also make it possible to set.- 
up larger loops by which implicit functions may be computed. Larger loops, how- 
ever, are not as ifoolproof and may become unstable, either oscillating or else 
changing as rapidly as possible to some value other than the desired one. 

The fact that the incremental computer is rate-limited becomes quite impor- 
tant when dealing with loops. This limiting may restrict changes in such a 
manner so as to diminish overshoot. On the other hand, there have been cases 
where it limited negative feedback so as to allow instability. In general, how- 
ever, the usual analytical methods developed for feedback loops should apply to 
incremental loops as well. 

Once a sequence has been selected, the next step is to determine incremental 
methods for generating it, ■ In. addition to the specified operations this genera- 
tion may involve such techniques as polynomial approximation and double integra- 
tion for sine and cosine. 

On long programs, especially, it is desirable to pack as much on to the drum 
as possible. This means performing as many operations in each step as possible. 
Though usually only one operation can be performed at a time, there are a few 
cases where two or more operations can be done at once. One of the most obvious 
cases of this is the use of the P-input, By using it, a term may be added :. * 
at the same time another operations is being performed. The chief limitation 
to the P-input is that it must be scaled the same as W e This usually is not a 
serious problem. Where it is, however, it may be necessary to use a separate 
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addition step. Luckily, addition steps usually have short word lengths so that 
they are not as costly space-wise as most other steps. In some cases addition 
can be accomplished using only the Q- and P-inputs, thus leaving the V-line 
available for comparisons. As a final remark it might be mentioned that repeti- 
tion should be kept to a minimum with no function being generated more than 
once if at all possible. It should be pointed out that it is less costly as a 
rule to add in a function already generated than to regenerate it, 

(2) SEQUENCING. - After the programmer has determined the set of func- 
tions to be used, the next step is to arrange them in the best possible sequence. 
In this operation the programmer aims '!) to have minimum lag between input and 
output, i.e., minimum computation time, and ;2) to efficiently utilize the RAM 
(random access memory) so as to have adequate memory positions available. 

In the present computer, due to /.I) the fact that the result of a computa- 
tion is not stored until the step following, and (2) the fact that the incre- 
ments must be read in the step preceding the computation, a result computed in 
cycle n cannot be used until cycle n + 2, In arranging the sequence of steps, 
to obtain the minimum computation lag, mathematically consecutive operations 
should therefore fall on alternate cycles. The (n ♦ l)th cycle would use the 
result of the previous major cycle, thus causing a major cycle delay. 

The simplest means of spacing the operations as described above is the 2- 
interlace. Witi the 2-interlace, consecutive operations are written in alter- 
nate drum positions with the first half of the program occupying, say, the even 
numbered positions and the second half occupying the odd-numbered positions, 
**• With this arrangement there is a two drum revolution delay between the beginning 
and end of the program, though all of the quantities are computed twice within 
this period, (In the first revolution the first half of the program is com- 
puted. In the second revolution the results of the first half are used as a 
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basis for the computations of the second half while the first half is being si- 
multaneously recomputed.) The effect of the two-interlace is to reduce the 
computation delay from one cycle per step to two cycles over-all. It has the 
advantage of being simple to use with a program written out in consecutive 
fashion. It is possible, liowever, to reduce the delay still further in many 
cases. It often happens that there are sequences of operations which can take 
place simultaneously. It is possible to interlace these sequences so that both' 
are performed in the same time it would take for one using the overall two- 
interlace mentioned above. It is often possible to divide the steps in a pro- 
gram into two categories? those rapidly varying and those slowly varying,, If 
the rapidly varying functions are placed on alternate steps, tand the slowly 
varying functions placed in between, the result is essentially a one-cycle time 
delay. 

To obtain the shortest possible time delay it is usually necessary to use 
a combination of the above methods,, because a method as simple as the two- 
interlace will reduce delay to two cycles (for a completely consecutive program), 
more elaborate techniques may often be unnecessary. In some cases it is possi- 
ble to further compensate for delay by extrapolating ahead a period equal to 
the delay. A one-cycle delay occurs when an input is operated upon during the 
same cycle in which it is taken in. This is due to the fact that the increments 
for an operation are drawn f^om the RAM in the step preceding and that the re- 
sults of a comparison are not available until the cycle following. To obtain 
the shortest possible input delay, the input should be first used during the 
second step after the comparison. 

Generally speaking, one RAM position is needed for each step and for each 
input. As, the number of cores is limited (in this computer to 64), RAM capacity 
becomes a limiting factor in long programs. To expand the capacity of the 



computer in this respect, a second head was installed on the R-line to allow 
the following command: 

Wj = + sgn R i _ 1 
This enables the programmer in many cases to obtain Aw without referencing the 
RAM. Without this command most memory positions would ; be always occupied , 
since sgn R formed by a step is needed as AW in the step on the next revolu- 
tion. With it the position usually need be occupied only until the last other ' 
step calling for sgn R, As in a typical program sgn R may be used only during 
a few steps immediately following the step in which it is generated, the same 
core may be used to store several values of sgn R in the course of a revolution; 
one after the other. In long programs, the programmer may have to sequence the 
steps so as to make maximum use of this feature. In most cases, however, it 
would appear the capacity of the RAM is adequate so that no special rearrange- 
ment is necessary. In a typical problem that uses the entire drum capacity, 
about half of the RAM is needed. 

d, SCALING. - The aim of the scaling operation is to generate the sequence 
of functions selected in the programming operation, with the minimum of error. 
There are two principal sources of error to be dealt with; J) round-off error 
and .2) overloading error. The former is due to the fact that a function cannot 
be expressed any more accurately, on the average, than to the nearest whole in- 
crement, while the latter is due to the error which is introduced when a func- 
tion changes more rapidly than one increment per cycle so that the step generat- 
ing it cannot keep up. 

** Unfortunately, these two errors are related in such a way that to reduce 

one usually increases the other. When a step is scaled up, which is the only 
way to decrease round-off error for the step, then the maximum rate of the out- 
put of the step is decreased so that overloading becomes more likely. Similarly, 
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scaling down a step to decrease the likelihood of overloading increases the 
round-dff error for the step. The two types of error are different in nature; 
whereas round-off error occurs all the time, overloading occurs only in certain 
cases, for example, during settling, or when a variable is changing at a high 
rate. Thus, it is a case of comparing an error which is always present with the 
probability of another error being present. In cases where scaling is critical, 
a compromise normally has to be made. 

In scaling programs where the problem is completely defined and the values 
of the functions throughout. the computation are always known, a straight analy- 
tical approach may be used, scaling the steps so that the rate of their output 
is never greater than one increment per cycle unless overloading is deliberately 
allowed to occur on occasion for the sake of reducing round-off error. 

In actual practice, however, it usually happens that the system is not com- 
pletely known; with only general information on the functions being available. 
Two examples of this would be: 1) where the inputs can vary over wide ranges 
and are different each time, and 2) where the system is so complex that involved 
analysis makes scaling exceedingly time-consuming. One solution to these prob- 
lems lies in simulation, or in trial computation, whereby a representative set 
of examples can be run with tentative scaling and the final scaling based on 
these. A second approach to these problems involves the use of approximations 
by which the analysis can be simplified. An example of such procedures is con- 
servative scaling. 

(1) CONSERVATIVE SCALING. - The ideal case of scaling is when all func- 
<* tions are always varying at the rate of one increment per cycle without any 
10 steps ever overloading. As the rate of a function in an actual case normally 
varies, a function usually varies at less than one increment per cycle on some 
occasions and at more than this rate (i.e., overloads the step) on others. 
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Though round-off error is less than one increment, overloading may cause an 
error of many increments. It might seem desirable, therefore, to scale a step 
so that it will never overload under any conditions. This is known as conserva 
tive scaling. It has the advantage of being relatively easy to do, and, pro- 
vided the round-off error is tolerable, yields satisfactory results. It has 
the additional advantage of providing the shortest possible settling time for a 
given input to a step. The assumption for conservative scaling is that the in- 
put to each step can vary at the maximum possible rate, i.e., one increment per 
cycle, and in the direction causing greatest output rate. For this method of 
scaling, the ranges of the functions need to be known. 

Conservative scaling usually applies to steps where the A P-input is not 
used. In all direct operations the assumption that P can vary at one increment 
per cycle automatically means that W must be able to vary at the same rate, 
without even considering the effect of the other inputs. The A P-input, there- 
fore, will be assumed zero in the analysis of the individual operations which 
follow,.* 

Addition 

SAWj iU AQi +V ATi 

s fe lu | * |v 



Multiplication 

SAW = UjAV^ ♦ V^jAU. 

° = | u |max + | v |max 
Integration 

«* SAW £ U.A T + U. ,AT 

^ • * i-l 

W S* 2lul 

x 5s d U Uax 

Division 

v i-l AU i £ °i^ v i + saw^ 
*S will be considered as positive unless otherwise stated. 
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V • £ u + s 

I |min ■»- | J max 

Isl < Ivl. - lul 

||l |min I | max 
Note that for division the minimum value of V is the limiting factor. It is 

possible to have cases where no S exists which will satisfy the equation. 

Square Root 

2U A U > SAW 

S < 2|U| . 
« I |min 

Exponential 

2UAT< SAW 

5 s ^| u |max 
Logarithm 

2UAT > SAW 

S < 2 1 U I . 

= I Imin 

In conservative scaling the assumption is made that all functions vary at 
their maximum rate (i.e., one increment per cycle). Though this may often be 
the case during the settling period, during other times functions are usually 
changing at some rate less than maximum, and often this rate has some definite 
limit below one increment per cycle. To give improved scaling, a higher order 
method has been evolved. This modification assigns to each function an e (effi 
ciency) factor defined as its maximum rate in increments per cycle. Scale fac- 
tors are then calculated as above with the exception that appropriate terms are 
multiplied by their e-factors. 

With this type of scaling the effect of settling time becomes more impor- 
"? tant. As was mentioned above, in conservative scaling a step will settle as 

vO 

x fast as the inputs to it. In a completely conservatively scaled program, for 

a, 

example, the program will have settled as soon as the inputs have settled so 
that the overall settling time is merely the maximum input settling time, (In 
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the present computer, with an input range of 1000 increments, the overall set- 
tling time would be 500 cycles maximum for a conservatively scaled program, 
assuming all of the initial values to be at midscale.) When it is assumed, 
however, that the maximum rates are less than one increment per cycle, then 
some steps are going to overload during settling when rates are one increment 
per cycle. Settling will then take an additional period of time after the in- 
puts have settled, during which the accumulated values on the R-lines in the 
various steps will be reduced to zero. 

At this point a brief discussion of settling will be introduced. The 
settling period for a step may be arbitrarily divided into two parts: input 
limited and output limited. Input limiting is said to occur when the inputs 
are changing in such a way as to cause the output to vary at less than one in- 
crement per cycle. Output limiting is said to occur when a step is overloaded. 
The total settling time for a step may be calculated by first determining over 
which periods it is input and output limited. The times for the two types of 
settling are calculated separately and then added together. As an example, the 
case of multiplication will be considered. The factors U and V will be assumed 
to start from zero and proceed at one increment per cycle to the positive values 
^max an( * ^max wn * cn are greater than S. Input limiting will occur for S/2 
cycles, at which time the product will be S/4 increments. Output limiting will 

then occur until the value UV/S is reached. The total settling time is there- 
fore . { S/2 + HV^^} = Ulj^i = E^aax + 3 cycles> In th . s ■ -particular 

case the delay caused by the input limiting was equal to one-half the period 
over which input limiting occurred. 

An output limited step tends to accumulate large values on the R-line, thus 
necessitating longer word length. In general, a step should be scaled up until 
it either takes too long to settle or else cannot keep up. If settling were 
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not a limitation, then all outputs would presumably be scaled to change at a 
maximum rate so that conservative scaling and its modification would be identi- 
cal for all steps following input steps. As it is, settling time often appears 
to be a limiting factor so that the modification is of value. 

Conservative scaling is inefficient in that to insure that no step ever 
overloads, the steps are usually scaled to run appreciably under full capacity. 
This situation can be improved by the simple expedient of scaling the steps up- 
ward to a value, say 30 percent, above the conservative value. 

Because of its advantage, i.e., that all steps always keep up, conservative 
scaling should be used in cases where the accuracy is adequate. It provides 
minimum settling time and minimum word length. Its use is essentially mandatory 
in steps involving integration, for if an integrand does not always correspond 
to the function it is being integrated with respect to, integration error will 
be introduced. In the incremental computer this error is cumulative, and remains 
until computation is re-initiated. This applies also to operations implicitly 
using integration, such as generation of the logarithm and the exponential. 
(2) WORD LENGTH CONSIDERATIONS. - When a program is so long that it 
necessitat' s use of the whole drum, the word length of each step must be taken 
into account, A step which is output-limited to any great extent accumulates 
large values on the R-line. This action is most likely to occur during the 
settling period ast after the computer is turned on. One of the advantages 
of conservative scaling is that output limiting never occurs, hence the number 
on the R-line is never greater than 2S. 

"f In cases when conservative scaling is not used some criterion other than 

in 

x the magnitude of S must be used to assign word length. One method is to simu- 

late a series of extreme cases as a basis for picking word length. A straight 

analytical approach may also be used. It should be emphasized that while 
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temporary overloading is usually permissible, overflowing of the R-line must 
n ever be allowed to occur; for, if the R-line ever overflows, permanent, almost 
invariably serious, errors are introduced, making any additional computation 
useless. 

(3) INPUTS - RANGE EXPANSION. - Inputs which vary over only a fraction 
of the range from zero to full scale may often be scaled up so that this region 
occupies an entire comparison interval. This is permissible because of the 
fact that the digital-to-analog converter converts only the first ten low order 
digits of a number, ignoring the others. The new regions are equal to the ori- 
ginal (12 to 1012), plus some multiple of 1024. 

(4) EXAMPLES SCALING x . - As an example consider the function x* f 
where -5 < x < + 3 and x < ,02/cycle. First x is scaled. Assume, as is the 
case with the present computer, that the input range is + 500 increments. As 
the maximum value of x is 5 it would appear possible to give x a coefficient 
of 100, thus, representing it over the range from -500 to +300 increments. To 
do this, however, would give the computer function a maximum rate of 1/C or .01/ 
cycle. This is one-half the maximum rate of the external function. Now, while 
this condition might be tolerable, e.g., when the peak rate occurred only rarely, 
in general it is desirable to have the computer function always keep up. This 
is accomplished by specifying C = 50. The range of the function in increments 
is now from -250 to +150, a total range of 400. If it is assumed that x may be 
at any value when computation is initiated, then the logical initial value for 
x would be -50. If more is known about the behavior of x, a more suitable 
choice may be made. 

^ Having scaled the input to the computer, the next step is to scale the func- 

tion x^ itself. First of all conservative scaling will be used. Referring to 
the algorism for squaring it will be noted that the quantity x appears in both 
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the U and V registers and that it is essentially added to R, every time x in- 
creases by one increment. The largest value which appears in these registers 
is -250, corresponding to x = -5. In the worst possible case, therefore, R 
changes by 500 every cycle, due to this action. To compensate for this we 

(C,,) 2 <fcn^2 
simply set S = 500. The coefficient of the output C w is — |— - = -^g^— = 5 

x 2 = (5) 2 = 25, 5(25) = 125 increments. Using a different form: 
max 

(-250) (-250) _ 19 „ 
(500) " L£a * 

In this case the same result may be obtained using the rate of the function, 

-5 = 2xx, ^ = 10(.02) = .2, 1/.2 = 5 
A A max 

As a further illustration assume that the coefficient of 5 is not large 
enough to provide the desired accuracy and that for that reason C was raised to 

p 

10, (S = 250). In this case the maximum value of Cx^ equals 250. If x varies 

x > 2.5, however, x* will not keep up. For all values 

p 

less than this (which in some cases may be most of the time) x will keep up. 
If it is known that x is small whenever x is large, then this type of scaling 
would keep up all of the time. 

Consider now the settling time of x^ under the two types of scaling, assum- 

o 
ing x to be at its maximum. For conservative scaling (C = 5), x will settle 

just as soon as x will; where x = - 1, in 200 cycles. In the second case, 

however, this is not true. Assuming that x Q = - 1 as above, x will keep up for 

125-50 = 75 cycles at which time x 2 = 6.25,, x will then change at its maximum 

rate (,1/cycle) until it reaches the final value of 25. The total settling time 

j§ is 75 + 10(25 - 6.25) = 263 cycles. 

^ e. COEFFICIENT RELATIONSHIPS. - This section deals with the algebraic equa- 

tions and methods for computing scale factors and coefficients for the various 
operations. 
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(1) CONVENTIONS. - Whenever a function f is represented in the computer 
it has two functions associated with it: the function F and the constant, or 
coefficient, C fc These three functions are related by the equation 

F = C f *f CD 

F is the computer representation of f and is defined as the number of increments 
representing f. The coefficient, Cf , relating the two has units of increments 
per unit. As an example, if a distance x is represented in the computer by 10 
increments per foot, then when x = 7 ft., X = 70 increments. C x = 10. It 
should be emphasized here that the basic algorism and all equations associated 
with it deal with the function F, not f. 

(2) BASIC ALGORISM AND ASSOCIATED EQUATIONS. - In the addition opera- 
tion the increments A U and A V are programmed to be zero so that U and V re- 
main at their initial values U and V Q and act as scale factors. It is not the 
actual values of U, V, and S that are important, but rather their ratios. This 

is fairly evident from the basic addition algorism 2. 

U Q + V T U v 
w s P = S^ Q * s^ T * P (2) 

To obtain the coefficient relationships begin with the general equation 3. 

w = kq + mt + p (3) 

By definition 

P = C pPf U = C u u, V = C v v, W = C w w, Q = C q q, and T = C t t (4) 

Substituting equation 4 in equation 3 

W _ 0_ T P_ 
u w u q u t L p 

T ^w ^w ^w 

>o orW = kc^O+mc;T + c^P (5) 

x 

From equation 2 and equation 5 

!!fi-*£s ■ ' V o C w 

S " C <T = m C7. and c w = c p (6) 

q t 
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r - D o £a - v » c t - c 
or c »-sk -s";r- C P 

or * ~ T~ r F~ 

k t w m l> w 

The usual procedure for picking U , V Q , and S is to pick the smallest inte- 
gers yielding the desired accuracy for the ratio. Small integers are picked to 
keep the remainder as low as possible. As for all direct operations, the P- 
term must have the same coefficient as W. 
Product 

Taking the multiplication algorism equation 7 and substituting equation 4 

into the general equation for multiplication we obtain 

UV 
W ± f- + P (7) 

w = kuv + p (8) 

L. _ k L. L. + P_ 

Cw C u Cy Cp 

k C 





W = w 

c u c v 


uv +2* P 

C P 


let 


C w = C p 

S - CuCv 
k C w 




or 


c w = c p = 


C u C v 
k S 


Quotient 






Proceeding as 


for multiplication 




u - S(W - 
V 


P) 




u ^ k(w - 


p) 



(9) 



(10) 
(11) 



b u L w L p L v 



let L>M 
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Square Root 

u = k V^w - p c u 2 = k 2 (w - p) (13) 

Substituting equation 4 in equation 13 

c u 2 ~ k K c p ) 



let 



^w ~ ^p 

U 2 , 2 /W-P 



Cu 2 = k VT7") CM) 

from the square root algorism 

U 2 = S(W - P) (15) 

from equation 14 and equation 15 

k 2 C u 2 (k C u ) 2 ^ V^C 
C w " C w 
Integration - Trapezoidal 



S = " = p— - , C u - — 1 - (16) 



W =| /udT (17) 

w = k J udt (18) 

Substituting equation 4 in equation 18 

w_ = k f 5. d ?- = -±_ f udT 



k c w 



or W = — ^ /udT (19) 

C u C t 



from equation 17 and equation 19 
«* 2 k C w 2C u c t „ ^uCt 



£ S C u C t ' " kC w '^~ kS 

cu Int eg ration - Point Slope 

W ■ = 1 / UdT 



S - v C. > C w - v c (20) 



(21) 
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U_ t 51 /C w k C t dW 
= k dt, 

k C.C. 



U = 



Vt dW 



C w dT 



From equations 27 and 29 

S _ fc C u C t _ 2k C u C t ^ 2k C u C t 
w u w 



Using equation 21 in place of equation 17 

S " C u C t ' b "kC w ll w" k S u ^ 

Integrati on - Reciprocal Integrand 

' ■ Q =f '/jf <23) 

q=k/^ (24) 

Substituting equation 4 in equation 24 

Q_ = k f dw/c w =1£m fdW 
c q J u7c7 c w J U 

or 0»^/«. (25) 

From equations 23 and 25 

S _ kC q C u Q _ 2k C q C u S^ 

2 c w ' C w • C q = 2k C u C26) 

Differentiation 

Using the formula without the stabilizing term 

u = k dH " ' (28) 

at 

Substituting equation 4 in equation 28 



C u " K dt/C t " C w dT i ^ V} 



2 = Cw ' S ~ C, • Cw ~ « (30) 



Exponential 

W = e 2T / S (31) 
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This is obtained from 

kw = f kudt where w » u = e* (32) 

Substituting equation 3 in equation 32 

r=frf I r J r/»« 11 < 33 > 

c w J C u Cx C u Cx J 
In order to have AW = AU, C w = C u , therefore W = U 

W = L- f UdT (34) 

C T J 

From equations 17 and 34 

§=£-, S = 2C t . C t =| (35) 

From the above it should be noted that C w is independent of C t , the only require 
ment for consistency being that C w = C u . The coefficient of w is determined 
from the relationship 

W = C w kw (36) 

Logarithm 

T = | log e W (37) 

The logarithm is generated as the inverse of the exponential, namely, from the 
equation 

w = J ud(kt) (38) 

where w = u = x and kt - log e x (39) 

substituting equation 4 in equation 38 

r = k/F-dr = F~r- / UdT (40) 

L w J L u L t L u L t J 
dividing by equation 17 

2 k c w 

But as C u = C w as with the exponential 

2 = L c t =0 S=^i (42) 

Once again C t is independent of C w or C u and is determined only by the scaling. 
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f e INITIAL VALUES 

(1) Introduction, - At the time when incremental computation is begun, 
all U and V registers must be set to a consistent set of values, which are 
called initial values. This section deals with the theory of computing these 
register values once the initial point has been determined. Re-examination of 
the basic algorism reveals that there is one term (usually SAW) which is af- 
fected by the previous sign of the remainder. At the beginning, however, there 
should be no previous sign of the remainder to consider; therefore, the initial 
values should be set so as to cancel out any effects of this term for the first 
cycle. 

The following helps to explain the effect of this term. Assume all inde- 
pendent increments are zero. The increment for the dependent variable (deter- 
mined on the previous cycle), is defined as plus one for the first cycle. This 
increment multiplied by some register quantity is added into the R register, 
causing a second increment to be generated opposite in sign to the first. This 
second increment on the next cycle will counteract the effects of the first 
increment, thus leaving the R-line at or near zero. All steps following this 
one will have a -1 for the result of this step, whereas it will be first a +1 
and then a -1 for all steps preceding it. In order to compensate for this dif- 
ference the registers of all steps following the step are set to an initial 
value one higher than normal. The initial output for these steps is calculated 
using this higher value. After the first cycle, compensation occurs and the 
output is equal to the unraised value. 

(2) COMPARISON. - V is updated prior to comparison. During the first 
cycle_A V = + 1 so that (V + 1) is compared. For this reason the V register 
should be set to a value one less than the desired initial value. 



If the register values in successive steps are each raised by an increment 
as specified anove, the values in the final step could be considerably different 
from the values based on the unraised functions. After the computer is turned 
on these steps will settle to the unraised values as expected, but only after 
several cycles, as the raised value is several increments different from the un- 
raised value. This difference can be minimized by choosing the sign of the func- 
tions generated in the various steps in such a way as to have successive offsets 

9 2 
compensate for each other. For example, instead of generating x, x^, and x + y, 

9 o 
one could generate x, -x , and x + y. In this case the second step would be 

2 2 
said to be flipped (from + x to - x^) . Flipping is most conveniently done by 

changing the signs of increments in such a way as to leave the register values 

unchanged. 

For arithmetic operations, flipping does not affect the final answer, but 
merely eliminates the initial delay described above In the case of integration 
and related operations, however, this initial delay can cause integration error 
which is permanent. For this reason flipping should be used in steps affecting 
an integral. In other cases the procedure is desirable though not necessarily 
worth the effort. 

g. METHODS 

(1) INTRODUCTION. - The theory of program preparation has been dealt 
with above. The following is a presentation of the techniques and conventions 
which have been used in the task of program preparation. While no claim is made 
that this system is the best one, it has worked satisfactorily and is the result 
of modification of several prior systems with which trouble was experienced 
over the past months. 

The results of the preparation routine are expressed in the form of three 
tabulations: (1) an algebraic listing of the dependent variables in each step 
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accompanied by a sign column, step number, scale factor, and coefficient; (2) a 
program tabulation giving the signs and addresses of each of the increments in 
each step, and (3) a constants tabulation listing the values of U , V Qfl and S 
for each step. Tabulation 1 is essentially for checking and is the basis for 
tabulations 2 and 3. 

The first step of programming is to determine a set of intermediate func- 
tions to be used in. computing the desired outputs. These are arranged in se- 
quence in accordance with the principles described in the section on sequencing 
and listed on Tabulation 1, the sign column being left blank. On Tabulation 1, 
the quantities appearing on the V-line are singly underlined. Those steps which 
have inputs or outputs have the V-line quantity doubly underlined. 

The scaling operation is then begun, using whichever method is desired. 
Work in general is done consecutively from step to step in sequence. By this 
time the initial point has usually been decided upon so that settling time can 
be taken into consideration. If the coefficients are being primarily determined 
rather than the scale factors, then the scale factors must be calculated, too, 
rounded off to the nearest integer, and used as a basis for calculating exact 
coefficients. In this way the round-off error for the scale factors can usually 
be eliminated. It often pays to go first through the program and set approxi- 
mate scale factors or coefficients, performing the exact calculation later. In 
the case of conservative scaling it is often desirable to calculate W fflax in a 
step and use this value as a basis for scaling succeeding steps. This is parti- 
cularly convenient in cases where the functions in many steps maximize under the 
same conditions. 

The next operation is to calculate the initial values. Presumably by this 
time the initial point has been determined from input and settling considerations. 
When flipping is to be done, it is desirable at this time to calculate the 
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complete set of scaled initial functions and to tabulate them. Then, when the 
initial values are calculated, the sign of the dependent variable is taken in 
such a way as to make the values of the U and V registers holding this function 
as close to the calculated value as possible. As an example, consider the func- 
tion x 2 where x = 5, C v = 100, and Co = *0. The input to step is Cx + 1 =501, 

A. X 

and the output is (501) = 251. The step is then flipped to give - x 2 so that 

v 200 

250 = 251 - 1 rather than 252 may be entered into the U or V register of the 
succeeding step using x 2 . Whenever a step is flipped a minus sign is entered 
in the sign column of Tabulation 1, Initial values are most conveniently calcu- 
lated using the scaled values of the functions as shown above. For, using the 
scaled values it is possible to go from one step to the next, using only the 
previous values along with the scale factors in the computations. If these 
scaled values are checked with the initial function values mentioned above, a 
double-check on the scaling operation is provided. For convenience, the flip- 
ping of a step is done in such a way as to leave the quantities in the U and V 
registers for the step unaltered. Flipping a direct operation consists of re- 
versing the sign of A P, AQ, and AT. Flipping an inverse operation consists 
of reversing the sign of AW and the dependent increment. 

(2) SIGN CONVENTION. - The dependent variable of a step is defined as 
the function listed on Tabulation 1, neglecting the sign column. It is this 
value which appears in subsequent U and V registers. The function sqn R multi- 
plied by the sign in the sign column is taken as a positive increment. Inverse 
operations, i.e., division, square root, differentiation, and logarithm, are 

^ considered to have a negative sign if A f = - sgn R, which is the convention 

i 
o 

w used in the section on theory. 

x 

e-. 

In setting up the convention it was attempted to have a system whereby flip- 
ping could be done with a minimum of after affects. Too, it was desired to have 
a rigid convention in which errors were readily apparent. 
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Again it should be stated that though various operations are stated sepa- 
rately, the programmer must take all into account simultaneously to do the best 
job of scaling. Usually the operations are repeated, at least to some extent, 
in an attempt to converge to the best possible program. 

h. CHEEKING OF COMPUTED PROGRAM 

(1) INTRODUCTION. - Tabulation 2 provides a convenient intermediate form 
for a program. This is primarily because a punched tape of the tabulation can 
be fed into the Univac Scientific Computer for simulation. This form also has 
the advantage that it is convenient for checking errors. Due to the standard 
form and convention adopted a great many of the common clerical errors can be 
either avoided or else quickly spotted. The checks which are applied are 1) 
address checks, and 2) sign checks. A careful check of Tabulation 2 can almost 
eliminate programming errors. In the checking operations involving either Tabu- 
lation 2 or 3, Tabulation 1. is used as the standard. 

The following are the checks which are applied to Tabulation 2. The incre- 
mental addresses are checked against Tabulation 1 to determine: 1) whether the 
right addresses were used, and 2) whether they were addressed to the appropriate 
increments. Knowledge of the normal layout of the operations (see Figure 2) 
makes the latter simple to perform. 

Due to the convention used, the signs of the register values are always for 
the positive function (that is, for +f, even though f may be negative). They 
are unaffected by flipping the step in which they occur. The signs of A U and 
Av may therefore be checked to see that they generate + U and + V, The signs 
of the other increments may then be checked to see that they generate a function 
with the same sign as specified in the sign column of Tabulation 1. 

As the U and V registers contain the operands as listed on Tabulation 1, 
they may be checked against the initial value tabulations. A fairly convenient 
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Operation 

Addition 

Multiplication 

Division 

Integration (t) 

Integration (p.s.) 

Input-Output 

Square 

Square Root 

Logarithm 

Exponential 

Differentiation 



U 










i 
in 



NOTE: Each color stands for one incremental address. The dependent incre 
merit is green. 



Figure 2. Tabulation No. 2 Check 
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checks on the initial values is obtained by calculating them first using the un- 
sealed inputs and the set of equations and then recalculating them just on the 
basis of the U, V, and S register values alone using the algorisms, 

(2) SIMULATION. - The above checks allow one to be reasonably certain 
as to the correctness of his program. The next step is usually simulation on 
the Univac Scientific Computer. This is convenient at this time because Tabula- 
tions 2 and 3 may be punched on paper tape and fed directly into the computer 
once the SIMIC (SIMulation, Incremental Computer) control tape has been run in„ 
While the SIMIC routine allows the use of only static points as inputs, it, 
nevertheless, has been very successful for troubleshooting purposes. 

The usual procedure is to pick a static point at which all of the inputs 
have changed. The point is run until all of the computer values are expected 
to have settled, with periodic dumps taken during the simulation as well as at 
the end. The values for the computer functions are calculated for this second 
point and compared with the values obtained with the dumps. The comparison 
shows up any program errors which have not been found earlier. If errors are 
found, it is possible to check for errors in steps following the erroneous one 
by determining whether or not the later values are consistent with the first 
error. Once the program has been checked out in this manner it is possible to 
run any other simulation which is desired, including one with dynamic inputs 
(by means of the DYSIMIC and POLYSIMIC routines). When the program has proved 
satisfactory it may be coded up for the incremental computer itself. 

i. SIMULATION 

(1) INTRODUCTION. - There are three simulation routines which may be 
used: SIMIC, POLYSIMIC, and DYSIMIC. The difference between the three involves 
(1) different incremental input procedures, and (2) different monitoring features, 
In the SIMIC program the inputs (equivalent to the analog inputs for the real 
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incremental computer) remain constant unless programmed to a different value,, 
While this limited type of input is very useful for eliminating program errors, 
it does not permit simulation under dynamic input conditions which normally occur 
in practice. In some cases the incremental program itself may be used to pro- 
vide the dynamic inputs. This is done by incrementally programming functions 
to provide the desired inputs. A POLYSIMIC program is a SIMIC program in which 
incrementally-computed polynomials in rime are used for inputs. The third possi- 
bility is DYS1MIC. This routine, independent of S1MIC 9 takes values for the in- 
puts, tabulated for regular intervals over the proposed simulation period and 
uses cubic interpolation to provide values to the computer every major cycle. 
Independent computer routines may be used to provide the input tabulations if 
desired. 

(2) SIMIC. - The SIMIC input program is composed of several sections, 
each introduced by a code word. These code words prepare the computer to act 
appropriately upon the data in the section. The sequence of these sections in 
an actual program is important only in that the data conditions for any desired 
simulation must be given before the COMPUTE command is given. The program is 
all in Flex-code. 

(3) PROGRAM. - The section is introduced by the word PROGRAM preceded 
by a shift up following a carriage return. It contains all of the incremental 
addresses from Tabulation 2. The information consists of one line for each step 
of the program, each line being composed of a 3-digit octal step number followed 
by six signed 3-digit octal addresses,, The computer ignores all material follow- 
ing the code word until it reads in a carriage return followed by 3 octal digits. 
During simulation the sign of the remainder of each step is considered to be 
stored in a register bearing the step number as an address. The addresses, in 
order from left to right, are: All f Av, AP t AQ, AT, and AW. Any increment 
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may be programmed zero by using "n" as an address. The sequence of steps with- 
in the program, constants, and input sections is immaterial; the step numbers 
determine sequence. The results of a comparison in a step are considered to be 
stored in a register with the same units and tens digit as the step, but with 
3 for the hundreds digit. There can be comparisons only up through step 077, 
There must be a line beginning with the number of the last step in the program 
plus one followed by the words "end of program". This causes the computer to 
begin the next major cycle. Finally, there must be a line beginning with the 
number 400 followed by six unsigned 3-digit octal step numbers. The presence 
of a step number causes the V-line of the step to be printed out. The addition 
of 300 (octal) to the step number causes the R-line of the step to be printed 
out. The program may contain up to 277 (octal) steps. A stop code (anywhere 
in the program section) signifies the end of the section and causes the computer 
to begin looking for the next code word. However, a carriage return at the end 
of each line (including the last) is necessary to cause that line to be stored. 

(4) CONSTANTS. - This section, introduced by the sequence; carriage 
return, shift up, CONSTANTS, contains all of the initial values and scale factors 
from Tabulation 3. There is one line for each step. Each line begins with a 
step number followed by 3 signed 6-digit decimal numbers; U Q , V and S, in 
that order. Each number, though integral, may be punctuated by a period between 
any digits. The section is ended by a stop code. Again, each line (including 
the last) must be followed by a carriage return. 

(5) OUTPUT HEADING. - This section is introduced by the sequences OUT- 
PUT, space, HEADING, carriage return. All material following the carriage re- 
turn is stored character by character. It is punched out before each section of 
results. A stop code following a carriage return signifies the end of this ma- 
terial. It is limited to 376 (octal) characters, including spaces, shifts up or 
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down, punctuation, etc. The output heading may be eliminated by programming an 
output heading with no characters in it. 

(6) INPUT. - This section is introduced by the sequence: INPUT. All 
material is ignored until a carriage return followed by a step number is read. 
Each line consists of a step number followed by signed, 6-digit decimal number 
with a period between the third and fourth digits. There must be an input line 
for each step containing a comparison. The section ends with a stop code follow- 
ing a carriage return. 

(7) CONTROL, - The above sections are used for loading information into 
the computer. The control of the simulation is essentially performed by the 
words "type (one space) spacing", "compute", "dump", and "restore". The word 
"type spacing", preceded by a carriage return, is followed by a signed, 6-digit, 
decimal number with a period between the third and fourth digits. During the 
simulation, this number gives the number of major cycles between print outs 
(actually punched). The word "cycles" may optionally be added after the number; 
it has no effect. The number, once set, remains the same until changed. The 
word "compute", occurs with the same format as "type spacing". The number in 
this case, however, specifies the number of major cycles to be simulated; it is 
normally a multiple of the "type spacing" number. The "compute" command must 

be repeated for each computation desired. 

A "dump" is a print-out of all the simulated registers in a program. It 
may be called for at any time by inserting the word "dump", preceded by a car- 
riage return and followed by a space. This instruction by itself will print 
out the contents of all the R-lines. If the space is followed by a "V" t all of 
the V- and R-lines will be printed out. If the space is followed by a "U", all 
of the U-, V-, and R-lines will be printed out. This print-out, as is the nor- 
mal print-out, is in flex-coded decimal. A special tape may be added to the 
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SIMIC master program which causes octal rather than decimal print-out. Th« 
code symbol " =s " will cause six inches of leader to be punched. This is often 
nsefol i®x identifying the various sections of a tajpe^ 

Mien it is desirei &® sJUnulate semml programs |n succession without re- 
loading the master tape; the "restore** txommand allows this to be done by .setting 
all of the R-lines to zero and by setting all of the increments to their initial 
value, which is plus one. Restoration is caused by the word "restore" preceded 
by a carriage return and followed by a carriage return and then a stop code. 

As mentioned above, stop codes signify the end of certain sections Cany 
section of variable length) and cause the computer to search for the next code 
word. Any uncalled for stop codes will cause the computer to stop and may be 
used for this purpose. The computer may be started again by pushing the START 
button. 

(8) COMPUTATION TIME. - As would be expected, the computation time is 
a function of: ,1) number of steps, 2) number of major cycles, and /.3). number 
of print-outs. As a rough guide, the computation rate may be taken as four 
cycles per second for a one-hundred octal-step program. Print-outs require ap- 
proximately one second per line. As this program operates almost exclusively 
off the drum, it is essential that an 8-interlace be used for most rapid compu- 
tation. However, there is little loss in speed when a 16-interlace is used. 
Use of 4-interlace will increase the computation time prohibitively. 

(9) DYSIMIC. - For convenience, DYSIMIC has been designed to use essen- 
tially the same program and constants sections as SIMIC, It is more restricted 
with respect to format, however. The computer recognizes the "PROGRA" of 
"PROGRAM", and then ignores all else until it read a a carriage return followed 
by three zeros. It then takes in the following information, ignoring step num- 
bers. It assumes the steps to be consecutive and increasing. This process 



74 



stops on the recognition of the "e" from "end of program". Search is then be- 
gun for the 400-line which may contain the addresses of up to six V addresses, 
i.e. r step numbers. The following carriage return triggers the search for the 
500-line into which the control section is condensed. The "500" is followed by 
four octal number in order* v l) total major cycles, 2) major cycles between 
input references [note: ;2) =1) for no references], ;3) major cycles between 
print-outs, and ;4) steps per major cycle [,4) » 76 (octal)J , 

The carriage return following the 500-line triggers the search for the 
"CONSTA" of "CONSTANTS". After a carriage return followed by three zeros the 
constants are stored. Once again the step numbers are ignored and the steps 
are assumed to be in consecutive, increasing order, U , V , and S must be pre- 
sent in that order. The numbers are signed, 6-digit and decimal. The period 
used in SIMIC is ignored and thus need not be present.. The "c" from "end of 
constants" is sensed and stops the computer, ready for the actual simulation. 

J.' DEMONSTRATION PROGRAM. - An example is included below as an illustration 
of programming. In it, the functions A A and AE are computed. (The A here 
does not indicate an increment.) In the following pages the equations, program, 
and sample calculations are given. In addition, a SIMIC input tape is repro- 
duced. 

This program uses double point-slope integration to generate sine and cosiss 
functions. The use of this' type of integration leaves the V-line free for in- 
put/output. The initial values of the inputs were chosen and the corresponding 
T and Rf calculated. The values used are: 
R = 5400 ft., -R s ♦ 400 ft./sec. , 
(A - it ) .=s E * E =s * ■'=* 
« *",.94 f - V 3 « 6 » W= 600 ft./sec., <V "- 1800) » 451 ft./sec. 
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for these: T = 2.3264 sec. 
R f a 4469.4 ft. 
AE = 7.3 mill {radians 



Ml 
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TABLE NO, 1 - EQUATIONS 
y a V Q - 10" 4 P R f CW cos A cos E ♦ 2030) 

f 2(R ♦ RT) 

v- 
I 

T = R f /V ' 

Q = !4T 2 

AA a JtX.- / C.06O« + ,13a ^WTslnA + W wAffDfE< 

cos E ¥ 4e4 x 10 4 



AE = ET ♦ cos E * <«<*<>• a + > 13 3> ^o w T sin E cos A _ W sin A 

R f V 5.58 x I0 4 
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TABLE NO. 2 - INPUTS 



Inout 


Description 


Unit 


Coefficient 


Minimum 


Maximum 


R 


sight range 


ft 


.16667 


o 


6000 


(A - IT) 


sight azimuth . 


rad 


500 


' -1 


♦1 


E 


sight elevation 


rad 


750 


2- 
~3 


2 
3 


o 

-R 


range rate 


ft/sec 


1.25 





800 


T* 


traverse irate 


rad/sec 


1750. 


-.3" 


♦.3 


E 


elevation rate 


rad/sec 


1750. 


-.3 


. -.a" 


V 


gun platform velocity 


ft/sec 


1. 





1000 


(V<> = 1800) 


muzzle velocity 


ft/sec 


1. 


1800 


2800 


*6 


relative air density 


«... 


1000 





1.0. 



Ca - .7833) relative air temperature — 



4615 



.5666 1, 



: j 



= A cos 
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TABLE NO. 3 - TABULATION NO. 1 
Step Sign Function Intuits 

R 



OOO 


♦ 


RT 




001 


♦ 


sin E cos A 




002 


♦ 


R ♦ RT 




003 


+ 


r2j2 




004 


♦ 


E 2 




005 


+ 


Q = 14T? 




006 l 


+ 


E)2 ♦ ±2 




007 


+ 


cos A cos E 




010 


4- 


(a - .7833) T 




Oil 


+ 


Q cos E/R f 




012 


* 


PT (.060 a + .133) 




013 


+ 


R 2 T 2 (f 2 + £2) 




014 


+ 


P TW (.060 a + .133) 




015 


+ 


R 2 T 2 (f 2 + e 2 )/2(R + RT) 




016 


* 


W cos A cos E 




017 


+ 


R f = R + R T + iMli^JLl?! 

1 2(R + RT) 




020 


+ 


Tl 




021 


+ 


sin E = f cos EdE 




022 


+ 


R f (W cos A cos E + 2030) 




023 


+ 


f T/cos E 




024 


■»■ 


ET + Q cos E/R f 




025 


•»• 


V = V - 10" 4 p R f (W cos A cos 


E + 2030) 


026 


+ 


cos E = - f sin EdE 




027 


•*■ 


P TW(060 a + _133)/V 





Note: Single underline indicates V-operand; double underline 
indicates input or output. 



-R 



W 



(V - 1800) 



(A - w) 
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TABLE NO e 3 - TABULATION NO. 1 - Ccont.) 
Step Sign Function Inputs 

R f / V a T 

P Bl sin A (.060 a ♦ .133)/ V ♦ f T/cos E = AA 
P TW sin E cos A (.060 a + .133)/ V ♦ Q cos E/R f ♦ 

ET = AE 

sin A a J cos AdA A A out 

(.060 a ♦ ,133)T * .06 [(a - .7833) T * 3t] 

cos A = - J 8 ln AdA AE out 

fi cos E 



030 


♦ 


031 


♦ 


032 


♦ 


033 


* 


034 


* 


035 


♦ 


036 


«»■ 



Note: Single underline indicates V -operand, double underline 
indicates input or output. 
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TABLE NO. 4 <* SAMPLE SCALING CALCULATIONS 
W cos A cos E I 

S- Wff 1 -' *800 

C = s8 ^5 Si = 1.00105 



rT 



s . C1750)X3 QQ1 5262 
5 " (2000) M 



C = 2003.8 



Rf/ V » T- 1 

C a (1800X1/6) = 300.00 
Rf (W cos A cos E * 2030) \ 

. S- (2000) ~ 83 - 4208x 10 
i(E) (T) «■ Q cos E/R f I 

c . U750)(30Q) , 20Q3a 
V (262) 4J003 * 8 

fT/cos E I 

« ■ S - <1214.90)<50Q) ~ 3Q3 

B (2003.8) * 303 

■ c g (3Q3) (2QQ3.fi) , 499 ?5 
h (1214,90) 4VV ' 75 
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TABLE NO. 5 - SAMPLE INITIAL VALUE CALCULATIONS 
»QQ RT | 

= 420 
*01 sin E cos A I 

W ~ (892) + l ' m * 

=* - 1 

*02 R + RT 



-"gg" 1 *""* 1 = 746 



= - 155 



*03 R 2 T 2 



(421) 2 
W = A ^- ♦ 1 = 267 

= 266 

*04 E 2 | 

(♦l) 2 
W * (875) * l = * l 

*05 Q = 14T 2 J 

153 



« (3214) 2 "" 

* a 152 
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TABLE NO. 6 - INCREMENTAL COMPUTER PROGRAM - TABULATION NO, 2 



Program No. 2032-B 

Title; Demonstration Program 



Date: 9/10/56 



Programmer: GAC 
|e JL o* JL 



Step 


A U 


A V 


A P 


A Q 


A T 


A W 


Remarks 


000 


+ 


030 


+ 


300 




n 


+ 


300 


+ 


030 


+ 


000 




001 


+ 


035 


+ 


021 




n 


+ 


021 


+ 


035 


+ 


001 




002 


+ 


030 


+ 


302 


+ 


300 


- 


302 


- 


030 


+ 


002 




003 


+ 


000 


+ 


000 




n 


+ 


000 


+ 


000 


+ 


003 




004 


+ 


324 


+ 


324 




n 


+ 


324 


+ 


324 


+ 


004 




005 


+ 


030 


+ 


030 




n 


+ 


030 


+ 


030 


+ 


005 




006 


* 


306 


♦ 


306 


+ 


004 


+ 


306 


+ 


306 


+ 


006 




007 


+ 


035 


+ 


026 




n 


+ 


026 


+ 


035 


+ 


007 




010 


+ 


030 


+ 


310 




n 


+ 


310 


+ 


030 


♦ 


010 




Oil 


+ 


017 


+ 


Oil 




n 


- 


Oil 


- 


017 


- 


036 




012 


+ 


034 


+ 


312 




n 


+ 


312 


♦ 


034 


+ 


012 




013 


+ 


006 


+ 


003 




n 


+ 


003 


♦ 


006 


+ 


013 




014 


+ 


012 


+ 


314 




n 


+ 


314 


+ 


012 


+ 


014 




015 


+ 


002 


"+ 


015 




n 


- 


015 


- 


002 


- 


013 




016 


♦ 


007 


+ 


314 




n 


+ 


314 


+ 


007 


+ 


016 




017 




n 


+ 


317 


+ 


002 


+ 


015 




n 


+ 


017 




020 


+ 


306 


+ 


030 




n 


+ 


030 


+ 


306 


+ 


020 




021 


+ 


026 


+ 


321 




n 


+ 


317 




n 


+ 


021 




022 


+ 


016 


4 


017 




n 


♦ 


017 


♦ 


016 


+ 


022 




023 


+ 


026 


+ 


023 




n 


- 


023 


- 


026 


- 


020 




024 


+ 


030 


+ 


324 


+ 


Oil 


+ 


324 


♦ 


030 


+ 


024 
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TABLE NO. 6 - INCREMENTAL COMPUTER PROGRAM- TABULATION NO. 2 (coat.) 



Program No. 2032-B Dates 9/10/56 

Title: Demonstration Program 



Programmers GAG 
Page JL. °* -JL 



Step 


A U 


A V 


L 


i P 


A Q 


A T 


A W 


Remarks 


025 


♦ 


312 


♦ 


022 


♦ 


321 


- 


022 


- 


312 


■»• 


025 


■ 


026 


+ 


021 


♦ 


326 




n 


- 


317 




n 


♦ 


026 




027 


♦ 


025 


+ 


027 




n 


'- 


027 


- 


025 




014 




030 


♦ 


030 


+ 


025 




n 


- 


025 


- 


030 


- 


017 




031 


♦ 


027 


+ 


033 


♦ 


023 


♦ 


033 


♦ 


027 


•f 


031 




032 


♦ 


027 


+ 


001 


♦ 


Oil 


♦ 


001 


+ 


027 


•f 


032 




033 


♦ 


035 


♦ 


031 




n 


" ♦ 


326 




n 


♦ 


033 




034 




n 




n 




n 


+ 


030 


+ 


010 


+ 


034 




035 


♦ 


033 


P 


032 




n 


- 


326 




n 


+ 


035 




036 


♦ 


026 


♦ 


005 




n 


.+ 


005 


♦ 


026 


4 


036 




037 




end of program 
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TABLE NO, 7 - INCREMENTAL COMPUTER PROGRAM - TABULATION NO e 3 



Program No, 2032-B Dates 9/10/56 

Title: Demonstration Program 



Programmer; GAC 



• JL of JL 



S 



Step 


% 


Vo 




s 


c 






000 


♦ 


000,699 


+ 


000.899 


+ 


001.500 


„ 03333 






001 


- 


000.594 


♦ 


000.001 


+ 


000.892 


500,48 






002 


4- 


000.699 


+ 


000.499 


+ 


002.250 


.16667 






003 


♦ 


000,421 


♦ 


000.421 


♦ 


000. 667 


1.6666 x 


io- 6 

8 


004 


♦ 


000.001 


♦ 


000.001 


♦ 


000.875 


3500 


fl 


005 


+ 


000,699 


♦ 


000.699 


+ 


003,214 


2.0002 




006 


♦ 


000, 001 


- 


000. 001 


♦ 


000.875 


3500 






007 


- 


000,594 


♦ 


001.216 


♦ 


000.903 


800.84 






010 


♦ 


000.699 


+ 


000.722 


+ 


001.153 


1200,78 






001 


+ 


000.746 


•♦• 


000.035 


♦ 


000.167 


2003.8 






012 


+ 


000.679 


+ 


000.599 


+ 


001. 100 


1399.5 


. 


. 


013 


+ 


000.001 


+ 


000.267 


+ 


000.250 


<£, «3«3«3fe X 


10~ 5 




014 


+ 


000.371 


+ 


000.799 


+ 


000,700 


2 






015 


+ 


000.746 


+ 


000.001 


+ 


000.596 


. 16688 






016 


- 


000.800 


* 


000.801 


+ 


000.800 


1.00105 






017 


+ 


000.001 


- 


000. 001 


+ 


000.001 


. 16667 






020 


+ 


000.001 


+ 


000.699 


+ 


000,262 


2003.8 






021 


+ 


001.215 


+ 


000.450 


+ 


001.215 


750. 00 




■ ' '■ . 


022 


+ 


001.230 


♦ 


000.751 


+ 


002.000 


83.421 x 


10" 6 




023 


+ 


001.216 


+ 


000.001 


+ 


000.303 


499.75 






024 


+ 


000.699 


+ 


000.001 


+ 


000.262 


2003.8 







85 



TABLE NO. 7 - INCREMENTAL COMPUTER PROGRAM - TABULATION NO. 3 (cont.) 



Program No. 2032-B Dates 9/10/56 

Titles Demonstration Program 



Programmers GAC 
2 of 2 



T 
2 



Step 


"o 


Vo 


s 


C 




025 


+ 


000.601 


■f 


000.460 


♦ 


000.834 


1.00024 




026 


+ 


000,001 


- 


000.001 


+ 


000.463 


1214.90 




027 


+ 


001.922 


♦ 


000.283 


+ 


001.500 


3000.0 




030 


+ 


000.699 


♦ 


001.922 


+ 


001.800 


300.00 




031 


+ 


000.283 


•f 


000.001 


♦ 


000.500 


3000.0 




032 


+ 


000.283 


+ 


000.001 


+ 


000.500 


3000.0 




033 


- 


000.595 


- 


000.001 


+ 


000.595 


500.00 




034 


•*• 


000.001 


+ 


000/012 


+ 


000,013 


1539.45 




035 


♦ 


000.001 


+ 


000.021 


+ 


000.420 


595.24 




036 


♦ 


001.216 


+ 


000. 153 


+ 


001.215 


2.0000 
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TABLE NO. 8 - DEMONSTRATION PROGRAM UNIVAC SIMIC INPUT TAPE 



Program No, 2Q32-B Date; 9/10/56 

Titles Demonstration Program 



rammers GTU GAG 
JL of m J Lm 



Step 


A U 


A V 


A P 


A 


A T 


A W 




000 


+ 


030 


# 


300 




n 


+ 


300 


* 


030 


+ 


000 




001 


+ 


035 


+ 


021 




n 


+ 


021 


+ 


035 


• + 


001 




002 


♦ 


030 


+ 


302 


+ 


300 




302 


- 


030 


♦ 


002 




003 


+ 


000 


+ 


000 




81. 


+ 


000 


+ 


000 


■ ♦ 


003 




004 


♦ 


324 


+ 


324 




n 


+ 


324 


+ 


324 


+ 


004 




005 


+ 


030 


+ 


030 




n 


+ 


030 


♦ 


030 


+ 


005 




006 


+ 


306 


♦ 


306 


■¥ 


004 


♦ 


306 


•¥ 


306 


♦ 


006 




007 


+ 


035 


- ♦ 


026 




n 


+ 


026 


* 


035 


♦ 


007 




010 


+ 


030 


♦ 


310 




n 


+ 


310 


<0> 


030 


♦ 


010 




Oil 


+ 


017 


+ 


on 




n 




on 


- 


017 


- 


036 




012 


+ 


034 


+ 


312 




n 


+ 


312 


+ 


034 


+ 


012 




013 


+ 


006 


+ 


003 




n 


+ 


003 


+ 


006 


+ 


013 




014 


+ 


012 


+ 


314 




n 


+ 


314 


•i- 


012 


+ 


014 




015 


+ 


002 


+ 


015 




n 


- 


015 


- 


002 


- 


013 




016 


+ 


007 


+ 


314 




n 


+ 


314 


+ 


007 


+ 


016 




017 




n 


+ 


317 


+ 


002 


+ 


015 




n 


+ 


017 




020 


+ 


306 


+ 


030 




n 


+ 


030 


+ 


306 


+ 


020 




021 


+ 


026 


+ 


321 




n 


+ 


317 




n 


+ 


021 




022 


+ 


016 


♦ 


017 




n 


+ 


017 


+ 


016 


+ 


022 




023 


+ 


026 


♦ 


023 




n 


- 


023 


- 


026 


- 


020 




024 


+ 


030 


+ 


324 


♦ 


on 


♦ 


324 


♦ 


030 


♦ 


024 
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TABLE NO. 8 - DEMONSTRATION PROGRAM UNIVAC SINIC INPUT TAPE Ccont.) 



j 

X 



Program No 2032-B Dates 9/10/56 

Title: Demonstration Program 



Page JL 



GTO GAC 
of JL 



Step 


A 


A V 


A P 


A Q 


A T 


A W 


■ 


025 


+ 


312 


•f 


022 


+ 


321 


- 


022 


- 


312 


♦ 


025 




026 


+ 


021 


♦ 


326 




n 


- 


317 




n 


+ 


026 




027 


♦ 


025 


♦ 


027 




n 


- 


027 


- 


025 


- 


014 




030 


♦ 


030 


+ 


025 




n 


- 


025 


- 


030 


- 


017 




031 


♦ 


027 


♦ 


033 


+ 


023 


+ 


033 


♦ 


027 


+ 


031 


• 


032 


♦ 


027 


+ 


001 


+ 


Oil 


+ 


001 


+ 


027 


+ 


032 




033 


+ 


035 


+ 


031 




n 


♦ 


326 




n 


-f 


033 




034 




n 




n 




n 


+ 


030 


♦• 


010 


•¥ 


034 




035 


+ 


033 


♦ 


032 




n 


- 


325 




n 


+ 


035 




036 


♦ 


026 


+ 


005 




n 


+ 


005 


+ 


026 


+ 


036 


V 


037 




end of prog i 


am 
















. 


400 




033 




035 




020 




030 




022 




036 




(STOP) 



























TABLE NO. 8 - DEMONSTRATION PROGRAM UNIVAC SIMIC INPUT TAPE (cont.) 



Program No. 2032-B Dates 9/10/56 

Titles Demonstration Program 



Programmers GTTJ GAC 
Page -J- of 5 



CONSTANTS 



Step 


"o 


Vo 


S 


-- -j 


000 


♦ 


000. 699 


♦ 


000.899 


+ 


001.500 




001 


- 


000.594 


+ 


000.001 


+ 


000,892 


I 


002 


+ 


000.699 


+ 


000.499 


+ 


002,250 




003 


+ 


000.421 


+ 


000.421 


+ 


000.667 




004 


4- 


ooo;ooi 


+ 


000,001 


+ 


000.875 


* 8 


005 


+ 


000.699 


♦ 


000.699 


+ 


003.214 




006 


+ 


000.001 


- 


0Q0.001 


+ 


000.875 




007 


- 


000.594 


+ 


001.216 


+ 


000.903 




010 


+ 


000.699 


+ 


000. 722 


♦ 


001.153 




on 


+ 


000.746 


+ 


000,035 


+ 


000. 167 


■ 


012 


+ 


000.679 


, ♦ 


000.599 


+ 


001.100 




013 


+ 


000.001 


+ 


000.267 


+ 


000.250 




014 


+ 


000.371 


+ 


000. 799 


+ 


000.700 




015 


+ 


000. 746 


+ 


000.001 


+ 


000.596 




016 


- 


000.800 


+ 


000.801 


+ 


000,800 




017 


+ 


000. 001 


- 


000.001 


+ 


000.001 




020 


+ 


000.001 


+ 


000. 699 


+ 


000.262 




021 


■f 


001.215 


+ 


000.450 


+ 


001.215 




022 


+ 


001.230 


•f 


000.751 


+ 


002.000 




023 


+ 


001.216 


■»• 


000.001 


+ 


000.303 


■ 


024 

_ 


+ 


000.699 


•*• 


000.001 


+ 


000.262 
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TABLE NO. 8 - DEMONSTRATION PROGRAM UNIVAC SIMIC INPUT TAPE (cont.) 



£ 



Program No. 2032-B Dates 9/10/56 

Title: Demonstration Program 



Programmers GTU GAG 
Page 4 of 5 



CONSTANTS 



Step 


Uo 


v 


- 
S 


- 


025 


+ 


000.601 


■+ 


000.460 


+ 


000.834 




026 


+ 


000.001 


- 


000.001 


+ 


000,463 




027 


♦ 


001.922 


•f 


000.283 


+ 


001.500 




030 


+ 


000.699 


•f 


001.922 


+ 


001.800 




031 


+ 


000.283 


•*• 


000.001 


+ 


000.500 


• 


032 


♦ 


000,283 


+ 


000.001 


+ 


000.500 




033 


- 


000.595 


- 


000.001 


+ 


000.595 




034 


+ 


000.001 


+ 


000.012 


+ 


000.013 




035 


•f 


000.001 


+ 


000.021 


+ 


000.420 




036 


+ 


001.216 


+ 


000. 153 


+ 


001.215 




(STOP) 
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TABLE NO. 8 - DEMONSTRATION PROGRAM UNIVAC SIMIC INPUT TAPE (cont*) 

Program No. 2032-B Dates 9/10/56 Programmers GTU GAC 

Title: Demonstration Program Page 5 of 5 

INPUT 

000 + 000.582 

002 + 000.645 

006 + 000.036 

010 + 000.723 

012 * 000.600 

014 * 000.800 

017 ♦ i000v327 / 

021 + 000.451 

024 - 000.092 

026 - 000.470 

CSTOP) 

OUTPUT HEADING 

RESULTS 

Program 2032-b Demonstration Program GTU GAC 9/10/56 

CYCLE dA dE T V Rf Q 

CSTOP) 

type spacing 000,050 

compute 000.800 cycles 



<a dump u 
S CSTOP) 
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3. CODING 

a. INTRODUCTION, - The solution to a mathematical problem is expressed as 
a series of equations. In order for a computer to solve a problem, it must be 
given instructions so that it will solve these equations. The translation from 
equations to machine instructions is called coding. 

All of the information necessary to code for the incremental computer is 
contained in the command code sheet (^Figure 3 ) and the algorisms. The follow- 
ing paragraphs will explain and illustrate their use, 

b, PROCEDURE IN CODING. - The equations are broken down into a series of 
simpler equations, called steps. Each step must be capable of being handled by 

one of the special algorisms. For example, the multiplication algorisms is 

UV 
W = — ♦ P, so two quantities could be multiplied and a third added to the pro- 
duct in one step. 

Using the special algorism, identify the variables of each step with the 

corresponding variable in the algorism. For example, let us compute the equa- 

UV 
tion Z = XY. From the multiplication algorism we are given W =-<r- + P. To 

change our equation to look like this we can write Z = XY + 0. It is now plain 
that P = 0, X = U, Y = V, Z = W. 

We have nqw identified the variables in our simple equation with the varia- 
bles in the algorism. Using the restrictions for the multiplication algorism we 
have 



a V av«i-%."h- {risiis:! 



ve 
ve 



So we instruct the machine to read + A and + AV from the address where AY 
is stored, and ♦ 4T and + AU from the address whe"re AX is stored. We can 
make P = by not reading it. 

When Rj is positive, the output increment of the step read by AW is posi- 
tive, and when Rj is negative the increment is negative. Therefore, if we 
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COMMAND CODES 



CO 
CM 

CO 
>+ 



p 3 


2 

2 

2 

2 

2 

2 

2 

2 

2 

2 - 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 



Figure 3. Command Codes, Incremental Computer 



P 2 


2 
3 
4 
5 
6 
7 
e 
e 
e 
e 
e 
e 
e 
d 
d 
d 
d 
d 
d 
d 


1 
1 







2 2 
22 



D. S. Address 

No Action 

Initiate Comparison Non-complement . . 

Initiate Comparison Complement .... 

End pulse. 

End pulse 

Initiate step ..... 

Initiate step 

Read + AP from last indicated address 

Read + AW from last indicated address 

Read + A V from last indicated address 

Read +AQ from last indicated address 

Read + A T from last indicated address 

Read + AD from last indicated address 

Read + A W from Sign of R^ ... . . 

Read - AP from last indicated address 

Read - A W from last indicated address 

Read - A V from last indicated address 

Read - AQ from last indicated address 

Read - AT from last indicated address 

Read - AU from last indicated address 

Initiate cycle . . . 

Comparator address 

1 

2 

3 
i 

t 

15 

16 Comparator address 
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COMMAND CODES 



Pn 










2 








2 


2 




2 


3 




2 


4 




2 


5 




2 


6 




2 


7 




2 


e 


1 


2 


e 


2 


2'- 


e 


3 


2 


e 


4 


2 


e 


5 


2 


e 


6 


2 


e 


7 


2 


d 


1 


2 


d 


2 


2 


d 


3 


2 


d 


4 


2 


d 


5 


2 


d 


6 


2 


d 


7 







? 







2 2 




1 


2 2 




1 
i 


2 2 




t 

7 
7 


2 2 



D. S. Address 

No Action 

Initiate Comparison Non-complement 

Initiate Comparison Complement 

End pulse 

End pulse 

Initiate step . . . 
Initiate step . . . 
Read + A P from last indicated address 
Read + AW from last indicated address 
Read + A V from last indicated address 
Read +AQ from last indicated address 
Read + A T from last indicated address 
Read + A U f 1 " ™ l ast indicated address 
Read + AW from Sign of R^.j ...... 

Read - AP from last indicated address 
Read - AW from last indicated address 
Read - A V from last indicated address 
Read - A Q from last indicated address 
Read - AT from last indicated address 
Read - AU from last indicated address 

Initiate cycle . . . . 

Comparator address ...... 

1 
2 
3 



COMMAND POSITION 
09 of x 
09 of x 

12th digit after Initiate Comparison 
12th digit after Initiate Comparison 
00 of x 

00 of x 

02 of x - 1 through last digit of x - 1 

01 of x - 1 through 00 of x 

02 of x - 1 through 00 of x 
04 of x - 1 through 02 of x 

03 of x - 1 through 00 of x 

03 of x - 1 through 01 of x 

p ^ 20 of x - 1 where p = sum of the number 

of digits in x - 1 and x. .p t 21. 
02 of x - 1 through last digit of x - 1 

01 of x - I through 00 of x 

02 of x - 1 through 00 of x 

04 of x - 1 through 02 of x 

03 of x - 1 through 00 of x 

03 of x - 1 through 01 of x 

digit period preceding 

00 of x 
digit* period following 
Initiate Comparison 



NOtES 

1. Command Position refers to the allowable 
digit positions for commands for the xth 
step. 



2. Address for storing 
sign of R. . . . 



05 of x + 1 



3. Address for storing results 

of the comparison . . . . 4th digit 

following end 
pulse 

4. e - even numbers j 
d - odd numbers 

5. Information from a digit storage address 
remains available until another D. S. 
address is programmed. 

6. Where command codes permit, combinations 
of commands may be programmed simultane- 
ously, and comparator addresses may be 
programmed with commands or O.S. addresses 
if their common parts can be made identical 

7. 2 C ~ 0, 1, 2 or 3 in this position 
2 2 = 4, 5 t 6, or 7 in this position 

8. DS - Digit Storage 

9. M.S.D. of the IT and Uj preceding initiate 
cycle must be 0-* 

10. First two M.S.D.'s of the Rj preceding 
initiate cycle must be 0. 



15 
16 



Comparator address 



digit period following 
Initiate Comparison 



I 

v© 

a 



program AW to be plus (+)AW, the required case is satisfied. To make AW = - 1 
when R| is positive, we would read - AW e To do this in the case of multiplica- 
tion would be incorrect, however, as the value for Rj would change away from 
rather than toward zero. 

The command code sheet gives a code number for every command and the allowed 
digit positions for that command to be given. For example, let us read one of 
the variables in the above equation; take + AV. The code number to read ♦ AV 
is given as 2e3, The e can be any even number through 6, The allowed posi- 
tions are "02 of x-1 through 00 of x", The meaning of this is made clear in 
note No, 1, The "x" refers to the step number with which we are concerned. The 
02 and 00 refer to digit positions within a step. The digit positions are num- 
bered starting with the initiate step command as zero; this is necessarily so, 
since the digit position for initiate step (IS) command is 00 of x. Thus, we 
see that +AV must be read in the step before the step in which it is to be used 
in computation, The same is true of the other variables with minor differences 
in initial digit positions allowed. 

After all of the variables have been read from their respective addresses, 
computation takes place automatically. All that remains is to store the sign 
of the remainder for use in the next major cycle. According to the sheet this 
is to be done at 05 of x * 1, the fifth digit position of the following step. 
Let us illustrate with a typical equation and how it is coded. 
Z = (A/d) log A - b 
1. Break into steps 

Step No. 1 A/d 

Step No. 2 log A 

Step No, 3 

Step No, 4 (Step No. 2) (Step No. 1) - b 



m 

x 

a. 



Each step can be generated by the use of only one special algorism. Note 

that step No. 3 is blank. This is because the next step after step No. 2 uses 

the results of step No. 2, The results of step No. 2 are not available until 

after they are stored in 05 of step No. 3; if read out as soon as possible in 

step No. 3 it will not be available for computation until step No. 4. Now let 

us look at the program for this equation and the step-by-step analysis of the 

program. This is a general format that ail programs follow. Let us assume 

that A A is found at address 10, Ab at Oil, Ad at 012. 

Digit 
Remarks Period 

Initiate cycle (IC) 

Initiate step (IS) No. (0) 

AA address ' » (1) 

Read - W from last indicated address (2) 

Ad address (3) . 

read +AV from last indicated address (4) 

read -AQ from last indicated address (5) 

A/d address (6) 

read +AU from last indicated address (7) 

read -AT from last indicated address ( (8) 

IS No. 1 (0) 

A address (1) 

+ Aw (2) 

+ AV (3) 

+ AU (4) 

no action (5) 

log A address (^) 

+ AQ (7) 

+ AT (8) 







Command 
Code 


Step No. 




2d7 
260 
010 
212 


Prepare Step No. 1 


-< 


012 
203 
214 
001 
206 
215 
260 
010 
202 


Compute Step No. 1 
Prepare Step No. 2 


< 


203 
206 
200 
002 
204 
. 205 
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Store Step No. 1 
Compute Step No. 2 
Prepare Step No, 3 



Store Step No. 2 
Prepare Step No. 4 



Compute Step No, 4 



« 

ft* 



Store Step No, 4 



< 



260 
200 
200 
200 
200 
001 
f 260 
004 
202 
Oil 
211 
002 
203 
204 
001 
206 
205 
260 
200 
200 
200 
200 
260 
200 
200 
200 
200 
004 



IS No, 2 

no action 

no action 

no action 

no action 

store sign Ri 

IS No, 3 

A/d [log A J - b addresi 

+ AW 

Ab 

-AP 

Store sign R2 (log A) 

+ AV 

+ AQ 

A/d 

+ AU 

♦■ AT 

IS No, 4 

No action 

No action 

No action 

No action 

IS No. 5 



(0) 
(1) 
(2) 
(3) 
(4) 
(5) 
(0) 
(1) 



(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

CIO) 

CO) 

CD 

C2) 

C3) 



CO) 
CI) 
C2) 
C3) 
C4) 
(5)' 



n 



I 



Explanation 

Step No. 2d7 



digit period 
00 



digit period 
01, 02 



Initiate cycle. This command begins operation. The d 

can be any odd number 1 through 7. "d M will hereafter 

be called 1. 

26 initiate step command goes at the beginning of 

every step. Any number can be put in the space. It 

will be called from here on. 

In the division algorism we have 
W - P 



U = 



VS 



•, To get our equation 



alsp 



*M = i 



P = 

AW = AA 

AV = AQ = Ad 

AU = AT = A(A/d) which is the answer or output of 

the step. 

+ 1 for Rj positive and Vj negative 
+ 1 for Rj negative and Vj positive 

- 1 for Rj positive and Vj positive 

- 1 for Rj negative and Vj negative 
This will be the result if -AW, - AQ, -AT are read. 

According to the command code we can get -AA on A W by giving first the 
address of AA(10) and then the command "Read -AW from the last indicated ad- 
dress" in digit period (dp) 01 or after. The address 10 was therefore put in 
dp 01 and the command "read -AW" in dp 02. 

dp 03, 04, 05 Similarly we read +AV and +AQ from address 12 where 

Ad is stored. Notice it is not necessary to repeat the 
address since the command reads "from last indicated ad- 
dress". 
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dp 06, 07, 08 ♦ AU and + AT are the output of the step, but must be 

read for use in computation. The place where the output 

of a step is stored is arbitrary, but the customary 

place is in the address whose number is the same as the 

step number in which it is computed; in this case 01. 

Step No. 1 dp 00 I.S. No. 1. Step 1 is automatically computed; step No. 

2 must be prepared. 

dp 01, 02, 03, From the algorism we have 
04 , 

Q == ^S log W 

AV = AU = AW = AA 

AT = AQ = Alog A 
A A is of course at address 10, and the following com- 
mands read it out. The address of log A is 002 as given 
by the convention mentioned in step 00. Digit period 
05 contains a no action command. If there were an ad- 
dress in this position the contents of the R line of 
step No. would be stored here and would erase anything 
previously stored. 
Step No. 2 dp 00 I.S, 

dp 01, 02, 03, There is nothing to prepare for step' No. 3, so the no 
04 

action command is given. 

dp 05 The sign of the remainder of a step is stored at 05 of 

x •*■ 1 according to the command code, so the remainder of 

*• . step No. 1 is stored at the address given here. 

^ Step No. 3 dp 00 I.S. 



nv 
dp 01, 02 To prepare for step No. 4, the algorism gives W = jj- ♦ P 

To form our equation 
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MA/d) = AU « AT, AP a -Ab, Aiog A a AV = 
The output of step No, 4 is stored at 04, so this is 
read by +AW in 01, 02. We want AP » - Ab, so these 
are read in 03 , 04. The sign of R2 is stored in dp 05. 
Log A is read from 02 by * AV and ♦AQ. Notice it is 
not necessary to repeat the address. A/d is read from 
01 where it was stored by +AU and + AT. 
Step No. 4 dp 00 I.S. There is nothing to be prepared or stored, and 

computation is automatic. The length of this step is 
somewhat arbitrary, but is set by factors which will be 
mentioned later. 
Step No, 5 The result of step No, 4 is stored in 05 of step No. 5. 
(1) CONSTANTS. - The U, V, and S constants are put into the computer in 
binary notation, A typical coding form might look like this; 



Constants 


Octal Program 
and Constants 


■ -■- 

Remarks 


~~5 


V Q 


u o 


P4 


P3 


P2 


PI 












2 


.6 





Initiate cycle 


1 






4 


- 


_ 


■_ 







1 




2 


- 


- 


- 




1 


1 





6 


- 


_ 


_ 
















_ 


_ 


_ 




1 





1 


5 


_ 


_ 


_ 







1 





2 


_ 


_ 


_ 




1 


1 


1 


7 
















1 


1 










1 








4 














1 


1 











I 
fiU 



The constants are staggered. This staggering causes S to reach the arith- 
metic section first, and forms SAP and -SAW. Next V arrives and VAT is formed 
and added to the rest, and so on. If we call the line containing initiate cycle, 
00? the least significant digit of S begins on 01, of V on 02, and of on 03, 
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E' 



Then the rows of digits are taken as octal numbers and the result is put In the 
P4 column, as illustrated above. Of course there must be enough digit positions 
in a step so that the constants can be put in. 

If a constant is a negative number, the 2°s complement of the positive con- 
stant is used. All complemented numbers have a series of Vs for the most sig- 
nificant digits. These l°s should be continued up to the constants of the follow, 
ing step. For example 

S V U P4 P3 P2 PI 



1 


1 1 


7 


1 


1 1 


7 


1 


1 1 


7 


1 


1 1 


7 


ol 


1 1 


3 


1 


~ol_L 


5 





i i 


3 


■1 





4 





1 


2 


1 


1 1 

etc. 


7 







Initiate step 



Thus, we see that each column is continued up to the following step, 
(2) INPUT/OUTPUT. - Comparisons^ 

Quantities are introduced into the incremental computer by means of compari- 
sons. In the process of comparing, the digital quantity in the computer on the 
V-line is converted to an analog current by means of a digital-to-analog con- 
verter. This current and the input current are run in opposition through wind- 
ings of a magnetic modulator. The output of the modulator is fed to a detector 
which generates an increment of either + 1 or - 1 depending on the sign of the 
current difference. In the case of computer outputs, the detector output is fed 
to a holding circuit which adjusts the value of the output current so that it 
approaches in magnitude the current from the digital-to-analog converter. 

The input circuitry is such that it can only accept currents in the range 
from zero to plus full scale. In order that variables which are both positive 
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and negative could be handled, a special command, complement comparison, was in- 
troduced. The inputs are biased externally so that zero for the variable corre- 
sponds to a half -scale current. Within the computer the highest order digit is 
complemented. The reason for this can be seen by considering the behavior of 
the highest order digit. Normally this digit is zero for numbers less than half- 
scale and one for numbers greater than half-scale. Using complement notation, 
however, the highest order digit(s) is zero for positive numbers and one for ne- 
gative numbers. If we adopt the convention that positive numbers be represented 
by currents greater than half scale, then it can be seen that the highest order 
digit must be complemented. 

The necessary commands for a comparison are: 

P3 P2 PI Command 

2 2 ■ ' - Initiate comparison non-complement 

- ' - - comparator address of modulator where variable appears 

2 '4 - end pulse 

- - address for storing results of comparison 

The following is a typical input step: 

Step (x - 1) 2 6 Initiate step (x-1) 

quantity appears on V-line 

2 6 Initiate step x 

1 
2 
3 
4 
5 
6 
7 

"? 8 

jjg 2 2 Initiate comparison 9 

- - comparator address 10 1 

E 2 

3 
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2 6 Initiate Step (x ♦ 1) 2 

11 
12 

2 4 end pulse -«•--•---• 13 

14 
15 
16 

» - address for storing 17 
results of comparison 

Explanation of Format . 

As previously mentioned, the quantity to be compared must appear on the va- 
line of the step during which it is compared,, Therefore, AV, is read during 
the step before,, i.e.* in x - 1, At 09 of x, the command "comparison non- 
complement" (or complement as the case may be) is given. The very next digit 
period, dp 10 of x, the address of the modulator which holds this variable Is 
given. The modulator addresses are given at the bottom of the command code 
sheet. For example, if our variable, A, were stored in modulator 3, we would 
give, in 10 of x, the command -14, -15, -16, or -17; they all have the same ef- 
fect here according to the command code sheet. Thirteen digit periods later the 
end pulse is given; this stops the comparison. Four digit periods later the 
address at which the result of the comparison is stored is given. This is the 
same as the address from which AV was read. In our case, A was stored at 010. 
This gives the A A needed in the computations and up-dates the value of A. 

Output , 

Output is exactly like input except the result of the comparison is not 
stored. As mentioned above, the result of the comparison is fed into a detec- 
tor and ultimately adjusts the output current. 

At this time it is thought that there should be a minimum of approximately 
35 digit periods between initiate comparison commands, 

(3) COMBINING COMMANDS. - According to the command sheet, commands may 
be programmed either simultaneously or with comparator addresses if their common 
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parts can be made identical. An example of this is a situation where the twelfth 
digit after initiate comparison came at the digit position where "read + AP" 
is. Both commands could be given by 241, 

(4) WORD LENGTH, - The number of digit positions in a step, also known 
as word length, must be at least as large as the largest of the following: 

1) The number of commands and addresses 

2) The number of digits on the R-line when it is at its maximum. 

If the R-line does not overload, its maximum is about 3S (assuming S>b, S>V); 
if it does overload e its value will be higher. If the word length is not long 
enough to accommodate the R-line, some of the most significant digits are lost 
and large errors introduced. 

(5) USE OF DOUBLE HEAD ON R-LINE. - Since there are 100 (octal) cores 
for storage of the output of a step or result of a comparison, the number of 
steps plus the number of comparisons would be limited to 100 octal (64 decimal) 
except for a special command. This special command allows the sign of a remain- 
der to be read from a special track, thus freeing a core for some other use. 
The command necessary for this operation is "Read + AW from sign of Ri.j", and 
the allowable digit position is P-20 of x - 1 where P = sum of digits in x - 1 
and x; P '£ 21, 

In normal practice this special command is used as often as possible if there 
are more than 64 steps plus comparisons so that core storage is available as 
needed during coding without re-arranging the command structure. 

When this special command is used practically the only limitation on the 

"* length of a program is the total number of digit positions. If the double head 

oft 

x on the R-line is used, the maximum word length is less than or equal to 20, and 

the sum of any two consecutive words must be greater than or equal to 21, 
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(6) STABILITY COMPARISONS. - In order to stabilize the modulators, pro- 
vision must be made in every program for zero and full scale comparisons. This 
is done in the following way. 

Zero Comparison . 

In some stepswhere the V-line is free, make the initial constant on the V- 
line equal to zero and do not program AV. This will keep V at zero. Then make 
an output comparison in the usual manner. At the time of this writing the modu- 
lator used for zero comparison is modulator 2, so the corresponding comparator's 
address would be put in the program. 

Full Scale Comparison . 

In some other step where the V-line is free, make the 10 least significant 
digits of the initial constant on the V-line equal 1. Then do not program AV 
and proceed as in zero comparison. The present full scale modulator is No. 1. 
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4. SIMULATION 

In order to test the behavior of the incremental computer, and to acquire 
facility in programming and scaling for the incremental computer, a program which 
simulates the logic of the incremental computer was written for the Remington 
Rand Univac Division's 1103 large scale computer. This section describes this 
program. The commands which comprise the program are listed at the end of this 
section. Reference is made to the following publications, obtainable from Rem^ 
ington Rand Univac Division, which describe the programming of the 1103. 

1) "Notes on the Logic of the ERA 1103 Computer" 

2) "The ERA 1103 Computer System", PX 71920, Vol. IV, Section 6, 
Programming. 

This simulation program was given the name SIMIC for SIMuiation Incremental 
Computer. It requires an input tape, the preparation and format of which have 
been previously described. Since this input tape contains constant inputs (R, 
A, E, etc.) SIMIC gives a static simulation of the Incremental Computer. A simu- 
lation program to receive a new set of inputs every major cycle is being pre- 
pared and has been given the name DYSIMIC e for DYnamic SIM uiation. Incremental 
Computer. 

The SIMIC program contains room for more computation steps (minor cycles) 
than the Incremental Computer itself. Consequently, a program to generate varia- 
ble inputs to the program may be included in SIMIC ahead of the program to be 
simulated. When used in this way, it is called POLYSIMIC since the variable in- 
puts are approximated by polynomials. This has been used primarily because DY- 
SIMIC was not yet ready. When it becomes ready, POLYSIMIC will no longer be 
used. 

When a simulation is to be made, a bioctal punched tape containing the SIMIC 
program is loaded into the 1103 computer by means of the Ferranti loading routine 
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in the Service Library. It goes on the drum in addresses 40000 through 44177. 
Drum addresses 44200 through 50000 plus 40x, where x is the number of minor 
cycles in the simulated program, are used for storage in the course of simula- 
tion The input tape is then inserted in the tape reader, PAK is set to 40000 
and the machine started. The simulation proceeds as per instructions on the in- 
put tape until the end of the input tape. If the code word "restore" appears 
at the start of a second input tape, this second simulation can be run without 
re-loading the SIMIC program tape. 

The SIMIC program can be divided into three tasks which the program must 
performs 

1) Read input tape; recognize code words; store constants and addresses 

2) Carry out the simulation 

3) Convert to decimal, flex-code, and punch out the results when called 
for by the input tape. 

The SIMIC program works from high-speed storage (HSS) in accomplishing tasks 1 
and 3 while 2 is done on the drum e 

Starting at 40000 the program jumps to 41776 which is the start of a block 
transfer. 40000 through 41777 is transferred to HSS and control passes to 00010 
(which now contains the instruction in 40010). This and the instructions follow- 
ing cause one character to be read from the input tape. This character is 
examined to see if it is a stop code (43) or carriage return (45). If it is a 
stop, the program jumps to 00007 which is an unconditional stop, with 40000 as 
NI. If a carriage return is found, the program jumps back to 00010 to read 

"f another character. If the character so examined is neither stop nor carriage 

in 

x return, it is stored in the scratch pad (00017 through 00037) and the program 

su 

jumps to 00100. As the characters are read in .they are stored in order in the 
scratch pad. Starting with 00100 the program examines the appropriate number 
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of characters in the scratch pad to determine the presence of one of the code 
words. If a given code word is not found, the program jumps to a section which 
seeks to identify the next code word. If a code word is found, the program jumps 
to the section which handles, in an appropriate fashion, the characters which 
follow that code word on the input tape. If none of the code words are found, 
the program returns to 00010, reads another character from the input tape, and 
again looks for one of the code words among the last n characters read, (n de- 
pends on the number of characters that comprise a code word). For example, the 
section beginning with 00100 seeks to identify the code word ''program". In this 
case n = 9 for the code word to be recognized includes "carriage return" and 
"shift up" characters as well as the seven letters of the word "program". The 
following table lists the code words, together with the first address of the 
section in which they are recognized, and where the program goes if the code 
word is or is not recognized. 



Address Code Word 

00100 cr-su-p-r-o-g-r-a-m 

00140 cr-su-c-o-n-s-t-a-n-t-s 

00200 o-u-t-p-u-t-sp, -h-e-a-d-i-n-g 

00240 t-y-p-e-sp. -s-p-a-c-i-n-g 

00300 -c-o-m-p-u-t-e 

00340 i-n-p-u-t 

01024 d-u-m-p-sp. 



If recognized If not recognized 



01223 
01252 



r-e-s-t~o-r-e 

where cr means carriage return code (45) 

su means shift up code (47)) 

sp. means space code (04) 

= means equal sign code (44) 



jump to; 


iunro tot 


42000 


00140 


00150 


00200 


01000 


00240 


01100 


00300 


01200 


00340 


00350 


01024 


01034 


01223 


01233 


01252 


01255 


00010 
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If none of the code words is found, control is returned to 00010 which reads in 
another character and initiates the above sequence again. This continues until 
a code word is found. 

When a code word is found, the program jumps to the indicated instruction 
where appropriate action is taken. When this action has been completed, the 
program jumps back to 40000 and, after the block transfer to HSS, begins to look 
for a "stop" or the next code word, A "stop" code at the end of the input tape 
(the incremental program being simulated) causes the computer to stop with 40000 
in PAK. 

The following paragraphs will sketch what happens when the various code words 
are recognized. The reader should have in mind the format of the input tape as 
well as the Incremental Computer algorism which this program is designed to 
simulate. 

When "program" is recognized, the SIMIC program jumps to 42000, This ini- 
tiates a block transfer of 42000 through 42777 into HSS, and jump to 00100, A 
character is read from the tape and checked to see if it is a relevant symbol, 
such as a number, carriage return, stop code, letters "e" (from "end of program") 
or "n", plus or minus sign, period, space, shift up, shift down, back space, 
delete, or tab. Any symbols other than these will cause the program to read 
the tape, ignoring everything, until a carriage return is found. When a number 
is found, it is converted from flex-code symbols into binary and the program 
jumps to 00140. Beginning at 00140 is a series of index jumps which cause the 
numbers to be organized into groups of three as they come in (see incremental 
.«*»'. tape format). There will be seven groups of three digits each, representing 

step, D, V, P, Q t T, and W, After each group is completed and stored, the pro- 

s 

gram goes to 00204. Here, and in the steps which follow, these addresses are 
manipulated In such a way as to produce the proper modifications in the steps 
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from 00400 to 00437. This block of steps is the nucleus of the whole SIMIC pro- 
gram, for it is here that the Incremental Computer's algorism is simulated. It 
may be indexed as follows: 

400 - 404 adds -SAWf to R^j 

405 - 411 adds Vj^ATj 

412 - 416 updates Us D A = V iX * A Uj 

417 - 422 adds ^ i t^Q i 

423 - 427 adds SAPj to form Rj 

430 - 434 updates V; Vi = Vj_j + AVj 

435 - 437 compares V. with input and stores the difference as the 
AV. to be used in this step at the next major cycle. 
After this block of instructions has been set up to carry out the algorism ac- 
cording to addresses given in, say, step 000, the carriage return at the end of 
the line signals that this task is completed. This causes this block to be 
transferred to the drum at 50000 through 50037 (for step 000) and then to get 
ready for the next line of addresses by again block transferring from 42000 
through 42777 into HSS. The next line is handled in the same way except that 
the algorism block (400 - 437) is transferred into 50040 - 50077 on the drum 
(for step 001). Thus, the algorism for each step is set up in accordance with 
the input tape and transferred to the drum. The instructions for step Xg appear 
on the drum beginning at (50000g ♦ 40qXq) , the subscript 8 denoting an octal 
number. When carrying out the simulation, the computation will start at 50000 
and run in order to (50000 ♦ 40y) where jump instruction has been placed by the 

*f program after finding the code word e-n-d (part of "end of program") opposite, 

to 

x say, step y. This is a jump to 00040. At the time the computation is carried 

out, HSS will contain instructions which decide after each major cycle, if cite is 

time to stopandi whether ;or;noti it lis < tlmettoppunchoout ^results. Thisccomes ifrom 
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the '"compute** and "type spacing" numbers on the input tape, and will be taken 
up later. 

After "end of program" , the input tape contains a line beginning with 400 
in the "step" column, (The six 3-digit octal numbers of this line, are the step 
numbers (of the incremental program) of those steps whose V-line values are to 
be punched out as output of the program.) When "400" is sensed, the SIMIC pro- 
gram jumps to 00500 where the instructions at 44005 through 44012 are changed 
to give the desired print-out. 

The "stop" code at the end of the "program" section returns control to 40000, 
and the machine begins to look for the next code word. 

Normally the next code word encountered is constants. The program jumps to 
00150, which now corresponds to 40150, Here, indexes at 43041 through 43043 
are restored and the program jumps to 43000. As before, there is a block trans- 
fer. The HSS now contains a duplicate of the instructions that are on the drum 
from 43000 through 43777. Starting again at 00100 the SIMIC program reads and 
decodes the step numbers and constants U , V Q , and S listed on the input tape, 
and stores these in their proper positions. These positions are listed in the 
table below. The table also includes, for convenient reference, the locations 
of other quantities used in the simulation. 

Ultimate Location 

in HSS at Time 

of Computation 

00200 - 00477 

00500 - 00777 

01000 - 01277 

01300 - 01377 
01400 - 01677 
01700 - 01777 



Quant it\ 




Location 


Stored 




on Drum 


V s 




44200 - 44477 


V 0i 's 




44500 - 44777 


V s 


< 


45000 - 45277 


Results 


of 




Comparisons 


45300 - 45377 


S 9 s 




45400 - 45677 


Inputs 




45700 - 45777 
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The location of a particular U 0f V Q| or S, within the range indicated above, de- 
pends of course upon the step number associated with it. Thus, SIMIC has room 
for up to 300 (octal) steps in the program to be simulated. SIMIC uses these 
locations as U-lines, V-lines, and R-lines are used in the Incremental Computer, 
i.e., the values in these registers are changed in the course of the computa- 
tion. To allow the initial values to be brought back, instructions located at 
43200 - 43224 mirror the U Q *s and V °s at 46000 through 46577, Later, if de- 
sired, the code word "restore" will clear the R-lines and return the U e s and 
V e s to 44200 through 44477, 

After the constants have been stored, stop code causes the program to go 
back to 40000 as before, and to look for the next code word. Normally this will 
be "Input", causing a jump to 00350 (40350). Here 43041 - 43043 are cleared 
and we jump to 43000. After the block transfer to HSS the program jumps to 
00100. Because the indices at 00041 - 00043 are now cleared, the program goes 
through the index jumps at 141, 142, and 143 (after each step number has been 
read in) and goes (via 00144) to 00210, Here, and in subsequent instructions, 
the input values are decoded, converted to binary, and stored in 45700 - .45777, 
depending on the step number. Again, the stop code at the end of the section 
causes a jump to 40000 and initiates the search for the next code word. 

The code word "output heading" causes a jump to 01000 (41000). The effect 
of these instructions is simply to store the flexowriter symbols, as read from 
the input tape, in locations starting with 41400, There is room for 376 (octal) 
symbols in the output heading, since any more than this would write over some 
of the instructions effecting a block transfer which is still needed. These 
instructions are located at 41775 - 41777, After storing the heading, control 
returns to 40000 and a search for the next code word is initiated. 
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The code word "type spacing" causes a jump to 01100 (41100) where the number 
following this code word is decoded, converted to binary, and stored in 44003. 
This time a stop code is not required to cause the program to return to 40000 
and seek another code word. 

The next code word is "compute". This causes a jump to 01200 (41200) where 
the number following the code word is stored in 44002. Then the stored output 
headings are punched out and the program jumps to 44000. 

The instructions starting at 44000 control the functions. First 44000 - 
45177 is block-transferred to HSS and control jumps to 00040. HSS now contains 
all the data needed by the computation which is on the drum starting with 50000, 
as well as the instructions which decide when to punch-out the outputs (type 
spacing) and when to stop the computation ("compute" number). Starting with 
00040, the program sets up the output registers, punches out the first line of 
outputs (for cycle 000), and jumps to 50000. Upon reaching 50000 + 40x (octal), 
where x is the step number corresponding to "end of program" on the input tape, 
one major cycle has been simulated. The jump to 00040 is made, and the cycle 
counter (at 00013) is compared with the "type spacing" and "compute" numbers. 
Whether or not an output is punched now, the program jumps to 50000 if the num- 
ber of cycles computed is less than the "compute" number. The program continues 
thus until these numbers are equal, when it again jumps to 40000 and looks for 
another code word. Usually this is a stop code, and the simulation stops. 

The remaining code words in the foregoing table are simply convenience fea- 
tures which have been added to SIMIC since its original inception. "Restore" 
has been described. Use of this word at the beginning of an input tape enables 
it to be fed into SIMIC after another simulation without the necessity of read- 
ing in the (large) SIMIC tape again. Used within an input program it enables 
the simulated program to • be brought back to its initial values of U Q and V and 
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run again with new input values. Thus, a dynamic simulation could be achieved 
by listing a new set of input values for each major cycle of computation. This 
is rather impractical, however, as it would mean a very long list of inputs. A 
"stop" code after "restore" is required to effect a return to 40000 and a search 
for the next code word. 

The "=" code word causes approximately six inches of leader to be punched 
out. This is convenient in identifying portions of the output tape. 

The code word "d-u-m-p-space" causes SIMIC to convert to decimal and flex- 
code, the contents of all the R*s in order. When immediately followed by "V", 
the program will punch out all the V's and all the R's, When followed by "U" 
(i.e., d-u-m-p-space-u) it will punch out the U's, V's, and R's. These dumps 
are invaluable in analyzing a proposed incremental computer program, especially 
in deciding upon scaling and word length which are governed in part by the 
build-up of the R-line. 

It is planned to add a "monitor" function to SIMIC also, but this has not 
been completely debugged as of the date of this report. Monitoring would com- 
pare the R-line of each step with three times that step's scale factor and 
punch out the cycle and the R-line when 35 is exceeded. It would do this each 
time an output is called for ("type spacing" number). The primary reason for 
adding this to SIMIC is that DYSIMIC, which has this useful feature, is not yet 
ready. It is planned to test Incremental Computer programs on POLYSIMIC until 
DYSIMIC is ready for use. 
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